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Abstract 

Gravitational waves (GW) are propagated perturbations of the space-time geometry, and they 
show up locally as tides, i.e., local gravity gradients which vary with time. These gradients are 
determined by six often called electric components of the Riemann tensor, which in a local coordinate 
frame are Romj, {hj = 1, 2, 3). A GW detector is a device designed to generate information on those 
quantities on the basis of suitable measurements. Both Weber bars and laser interferometers are 
single mode antennas, i.e., they generate a single readout which is the result of the combined action 
of those six GW's Riemann tensor amplitudes on the system sensor's output port. 

A spherical detector is not limited in that fashion: this is because its degenerate oscillation 
eigenmodes are uniquely matched to the structure of the above Riemann tensor components. This 
means that a solid spherical body is a natural multimode GW detector, i.e., it is fully capable of 
delivering six channel outputs, precise combinations of which completely deconvolve the six GW 
amplitudes. 

The present article is concerned with the theoretical reasons of the remarkable properties of a 
spherical GW detector. The analysis proceeds from first principles and is based on essentially no 
ad hoc hypotheses. The mathematical beauty of the theory is outstanding, and abundant detail is 
given for a thorough understanding of the fundamental facts and ideas. Experimental evidence of 
the accuracy of the model is also provided, where possible. 
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Summary 



The following pages contain the fusion of two published papers: What can we learn about gravitational 
wave physics with an elastic spherical antenna?, published in the Physical Review (PRD 52, 591-604 
(1995)), and Multiple mode gravitational wave detection with a spherical antenna, scheduled to appear 
in the July-2000 issue of Monthly Notices of the Royal Astronomical Society (MNRAS 316, 173-194 
(2000) )[]. In them, I first developed an analytic model to describe the interaction of a solid elastic 
sphere with an incoming flux of gravitational waves of the general metric class, then extended it to also 
address the multimode signal deconvolution problem by means of suitable layouts of resonant motion 
sensors. 

I think the model is quite complete as regards the major conceptual issues in a spherical GW 
detector, therefore think merging the two papers together makes positive sense. Let me however stress 
that I do not mean to understate the fundamental problem of noise, which is here left aside; rather, it 
is now being actively investigated in detail within the very convenient general framework set up below. 

What can we learn about GW Physics ... is already five years old, and certain parts of it have been 
the subject of further research since. More specifically, this happens with Brans-Dicke absorption cross 
sections, see section 4, where reference has been added to such newer work. In addition, two new tables 
and more mathematical detail have been added to its Appendix B. The latter should help the interested 
reader with certain technicalities, while the tables and attendant new formulas are meant to provide 
explicit numerical values of the frequency spectrum of the sphere rather than the original only graphics 
information. I expect this to be useful reference in numerical and/or experimental determinations of 
these quantities. 

No changes have been included in the second paper relative to its MNRAS version. In merging the 
two articles together, however, I have found expedient to renumber the equations and bibliographic 
references into corresponding single streams. The bibliography has been updated and made into a 
unique list at the end of the file in alphabetical order. Finally, I have included minor notation changes 
in order to reconcile otherwise small mistunings between both articles. 



Alberto Lobo 
Barcelona, June 2000 



1 The research contained in this second paper was actually complete by early 1996, and its main results first presented 
in a Winter School at Warsaw (Poland) in march 1996 [B9J, then in a plenary lecture at the Second Edoardo Amaldi 
Conference on Gravitational Waves at CERN in July 1997 J4lJ] . Essentially in its present form — see acknowledgements 
in page ^i] below — , the paper was submitted to PRD in 1997. It was rejected after the Editor, Dennis Nordstrom, had 
already accepted it, on the basis of a positive report. The reason for this opinion switch was a later, extremely vague yet 
demolishingly negative and utterly arrogant judgement by a second referee, to whom Nordstrom grants unlimited credit. 

An unchanged version of the manuscript was, again, positively reported on for MNRAS, and so the article will at long 
last see the light, three years too late. 
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Abstract 

A general formalism is set up to analyse the response of an arbitrary solid elastic body to an 
arbitrary metric Gravitational Wave perturbation, which fully displays the details of the interaction 
wave-antenna. The formalism is applied to the spherical detector, whose sensitivity parameters are 
thereby scrutinised. A multimode transfer function is defined to study the amplitude sensitivity, and 
absorption cross sections are calculated for a general metric theory of GW physics. Their scaling 
properties are shown to be independent of the underlying theory, with interesting consequences for 
future detector design. The GW incidence direction deconvolution problem is also discussed, always 
within the context of a general metric theory of the gravitational field. 



1 Introduction 

Spherical antennae are considered by many to be the natural next step in the development of resonant 
GW detectors. The reasons for this new trend essentially derive from the improved sensitivity of a 
sphere — which can be nearly an order of magnitude better than a cylinder having the same resonance 



frequency, see below and [ 13 1 — , and from its multimode capabilities, first recognised by Forward [22] 



and further elaborated in pi], 3C]. 



Although some of the most relevant aspects of detector sensitivity have already received attention 
in the literature, it seems to me that a sufficiently general and flexible analysis of the interaction 
between GW and detector has not been satisfactorily developed to date. This theoretical shortage has 
a number of practical negative consequences, too. Traditional analysis, to mention but an example, is 
almost invariably restricted to General Relativity or scalar-tensor theories of gravity; while it may be 
argued that this is already very general, any such argument is, as a matter of fact, understating the 
potentialities actually offered by a spherical GW antenna to help decide for or against any one specific 
theory of the gravitational field on the basis of experimental observation. 

I thus propose to develop in this paper a full fledged mathematical formalism which will enable 
analysis of the antenna's response to a completely general GW, i.e., making no a priori assumptions 
about which is the correct theory underlying GW physics (other than, indeed, that it is a metric 
theory), and also making no assumptions about detector shape, structure or boundary conditions. 
Considering things in such generality is not only "theoretically nice" — it also brings about new results 
and a better understanding of older ones. For example, it will be proved that the sphere is the most 
efficient GW elastic detector shape, and that higher mode absorption cross sections scale independently 
of GW physics. I will also discuss the direction of incidence deconvolution problem in the context of a 
general metric theory of gravity. 

The paper is organised as follows: section 2 is devoted to the development of the general mathe- 
matical framework, leading to a formula in which an elastic solid's response is related to the action of 
an arbitrary metric GW impinging on it. In section 3 the general equations are applied to the homo- 
geneous spherical body, and a discussion of the deconvolution problem is presented as well. Section 4 
contains the description of the sphere's sensitivity parameters, specifically leading to the concept of 
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multi-mode, or vector, transfer function, and to an analysis of the absorption cross section presented 
by this detector to a passing by GW. Conclusions and prospects are summarised in section 5, and two 
appendices are added which include mathematical derivations. 



2 General mathematical framework 

In the mathematical model, I shall be assuming that the antenna is a solid elastic body which responds 
to GW perturbations according to the equations of classical non-relativistic linear Elasticity Theory p2[ . 
This is fully justified because GW-induced displacements will be very small indeed, and the speed of 
such displacements much smaller than that of light for any foreseeable frequencies. Although our 
primary interest is a spherical antenna, the considerations which follow in the remainder of this section 
have general validity for arbitrarily shaped, isotropic elastic solids. 

Let u(x, t) be the displacement vector of the infinitesimal mass element sitting at point x relative 
to the solid's centre of mass in its unperturbed state, whose density distribution in that state is q. Let 
A and fi be the material's elastic Lame coefficients. If a volume force density f(x, t) acts on such solid, 



the displacement field u(x, t) is the solution to the system of partial differential equations [32| 



^-/xV 2 u-(A + /i)V(V-u) = f(x,i) W 

with the appropriate initial and boundary conditions. A summary of notation and general results 
regarding the solution to that system is briefly outlined in the ensuing subsection, as they are necessary 
for the subsequent developments in this paper, and also in future work — e.g. page ^ and ss. below. 

2.1 Separable driving force 

For reasons which will become clear later on, we shall only be interested in driving forces of the separable 
type 

f(x,i) = f(x) 5 (i) (2) 

or, indeed, linear combinations thereof. The solution to (|l|) does not require us to specify the precise 
boundary conditions on u(x, t) at this stage, but we need to set the initial conditions. We adopt the 
following: 

u(x,0) = u(x,0) = (3) 

where '= d/dt, implying that the antenna is at complete rest before observation begins at i=0. The 
structure of the force field @ is such that the displacements u(x, t) can be expressed by means of a 
Green function integral of the form 

/•oo 

u(x,*) = / S{-x;t-t')g{t')dt' (4) 
Jo 

The deductive procedure whereby S(x; t — t') is calculated can be found in many standard textbooks 



see e.g. [30|. The result is 



where 



if t < 

(5) 

J2n ift u N(^)sinuj N t if t > 



fN ^h I ^W-f(x)d 3 x (6) 

M J Solid 
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and Ujv(x) are the normalised eigen-solutions to 



/A7 un + (A + n) V(V-uat) = -oj n qun (7) 

with suitable boundary conditions. Here N represents an index, or set of indices, labelling the eigenmode 
of frequency un- The normalisation condition is (arbitrarily) chosen so that 

/ u%,(x) -ujv(x) g(x)d 3 x = M5 N , N (8) 

./Solid 

where M is the total mass of the solid, and the asterisk denotes complex conjugation. Replacing now (|5|) 
into (D we can write the solution to our problem as a series expansion: 

u(x,t) = J2—"N(x)9N(t) (9) 

where 

g N {t)= f g{t')smoj N {t-t')dt' (10) 
Jo 

Equation (||) is the formal solution to our problem; it has the standard form of an orthogonal 
expansion and is valid for any solid driven by a separable force like (^) and any boundary conditions. 
It is therefore completely general, given that type of force. 

Before we go on, it is perhaps interesting to quote a simple but useful example. It is the case of 
a solid hit by a hammer blow, i.e., receiving a sudden stroke at a point on its surface. Exam of the 
response of a GW antenna to such perturbation is being used for correct tuning and monitoring of the 



device [29]. If the driving force density is represented by the simple model 

f( hb )(x,t)=f ^ 3 )(x-x )^) (11) 

where xo is the surface point hit, and fo is a constant vector, then the system's response is imme- 
diately seen to be 

u (hb) (x,t) = V ^-u N (x)smu N t (12) 

N UN 

with f N = M~ l fo -u^(xo). A hammer blow thus excites all the solid's normal modes, except those 
perpendicular to fo, with amplitudes which are inversely proportional to the mode's frequency. This is 
seen to be a rather general result in the theory of sound waves in isotropic elastic solids. 

2.2 The GW tidal forces 

An incoming GW manifests itself as a tidal force density; in the long wavelength linear approxima- 



tion [54 1 it only depends on the "electric" components of the Riemann tensor: 



/i(x, t) = gc 2 Roi j(t) Xj (13) 



where c is the speed of light, and sum over the repeated index j is understood. In ( |13D tidal forces 
are referred to the antenna's centre of mass, and thus x is a vector originating there. Note that I have 
omitted any dependence of Roioj on spatial coordinates, since it only needs to be evaluated at the solid's 
centre. The Riemann tensor is only required to first order at this stage 



RoiOj — 2 (Aij.oo — h 0i j — hojfii + /ioo,ij) (14) 



5 



where h^ v are the perturbations to flat geometry^, always at the centre of mass of the detector. 

The form ( |i~3| ) is seen to be a sum of three terms like (||) — but this three term "straightforward" 
splitting is not the most convenient, due to lack of invar iance and symmetry. A better choice is now 
outlined. 

An arbitrary symmetric tensor Sij admits the following decomposition: 

S ij (t)=S^(t)E^ + S^ m \t)E^ m) (15) 

m=-2 

where E^ m ^ are 5 linearly independent symmetric and traceless tensors, and E®p is a multiple of the 
unit tensor dij. S^ lm \t) are uniquely defined functions, whose explicit form depends on the particular 
representation of the i?-matrices chosen. A convenient one is the following: 



E 
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The excellence of this representation stems from its ability to display the spin features of the driving 
terms in Such features are characterised by the relationships 

£g m > mrij = Y lm (n) , 1 = 0,2; m = -l,...,l (17) 

where n = x/|x| is the radial unit vector, and Y/ m (n) are spherical harmonics pi]]. Details about the 
above .E-matrices are given in Appendix A. In particular, the orthogonality relations ( |55|) can be used 
to invert fll5|): 




S m (t) = fE^S tJ (t) (18.a) 

S (2m) (t) = %E*^ Sij (t), m = -2,...,2 (18.b) 

where an asterisk denotes complex conjug ation. Note that S(°°)(t) = y/4nS(t)/3, where S(t) = f^S^) 
is the tensor's trace. 

We now take advantage of (|l^) to express the GW tidal force (|i~3| ) as a sum of split terms like (||) : 
f(x,t)=f(°°)(x)< 7 ( 00 )(t) + £ i (2m H*)g( 2m \t) (19) 

m=-2 

with 



2 Throughout this paper, Greek indices (fi, u, . . .) will run through space-time values 0,1,2,3; Latin indices (i, j, . . .) will 
run through space values 1,2,3 only. 
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/fV) = .4 0) *i . 2 (00) W = f 4 ( ° 0) W) c 2 (20.a) 

ft i \n) = e ^S m) x j , g^ m Ht) = ^Et^R mj {t)J (m = -2,...,2) (20.b) 

Straightforward application of (||) yields the formal solution of the antenna response to a GW 
perturbation: 



TV 



(21) 



with the notation of (|60 and ( |10|) applied mutatis mutandi to the terms in (|2C|). 

Equation (f2l|) gives the response of an arbitrary elastic solid to an incoming weak GW, independently 
of the underlying gravity theory, be it General Relativity (GR) or indeed any other metric theory of 
the gravitational interaction. It is also valid for any antenna shape and any boundary conditions, thus 
giving the formalism, in particular, the capability of being used to study the response of a detector which 
is suspended by means of a mechanical device in the laboratory site — a situation of much practical 
importance. It is therefore very general. 

Equation ( |2l| ) also tells us that that only monopole and quadrupole detector modes can possibly be 
excited by a metric GW. The nice thing about ([H]) is that it fully displays the monopole-quadrupole 
structure of the solution to the fundamental differential equations (|l]) . 

In a non-symmetric body, all (or nearly all) the modes have monopole and quadrupole moments, 



and (21) precisely shows how much each of them contributes to the detector's response. A homogeneous 
spherical antenna, which is very symmetric, has a set of vibrational eigenmodes which are particularly 
well matched to the form (ElJ): it only possesses one series of monopole modes and one (five- fold 
degenerate) series of quadrupole modes — see next section and Appendix B for details. The existence 
of so few modes which couple to GWs means that all the absorbed incoming radiation energy will be 
distributed amongst those few modes only, thereby making the sphere the most efficient detector, even 
from the sensitivity point of view. The higher energy cross section per unit mass reported for spheres 



on the basis of GR [13], for example, finds here its qualitative explanation. The generality of (|2l|), on 
the other hand, means that this excellence of the spherical detector is there independently of which is 
the correct GW theory. 

Before going further, let me mention another potentially useful application of the formalism so far. 
Cylindrical antennas, for instance, are usually studied in the thin rod approximation; although this 
is generally quite satisfactory, equation ( |2l| ) offers the possibility of eventually considering corrections 
to such simplifying hypothesis by use of more realistic eigenfunctions, such as those given in [27, 56]. 



Recent new proposals for stumpy cylinder arrays |8[ may well benefit from the above approach, too. 



3 The spherical antenna 



To explore the consequences of ( |2l|) in a particular case, the mode amplitudes u^v(x) and frequencies 
ujn must be specified. From now on I will focus on a homogeneous sphere whose surface is free of 
tractions and/or tensions; the latter happens to be quite a good approximation, even if the sphere is 



suspended in the static gravitational field [34] 



The normal modes of the free sphere fall into two families: so called toroidal — where the sphere 
only undergoes twisting which keep its shape unchanged throughout the volume — and spheroidal |H| • 
where radial as well as tangential displacements take place. I use the notation 

ul lm ^)e^ , u^x) e**** (22) 
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for them, respectively; note that the index N of the previous section is a multiple index {nlm} for each 
family; I and m are the usual multipole indices, and n numbers from 1 to oo each of the /-pole modes. 
The frequencies happen to be independent of m, and so every one mode (^) is (2/+l)-fold degenerate. 
Further details about these eigenmodes are given in Appendix B. 

In order to see what ( |2l| ) looks like in this case, integrals of the form (|6|) ought to be evaluated. It is 
straightforward to prove that they all vanish for the toroidal modes, the spheroidal modes contributing 
the only non-vanishing terms; after some algebra one finds 



fZm^Jjt I uP l * m ^)-i (l ' m ' ) ^)d 3 x = a nl 5 l/l 5 mlm , l' = 0,2, m' = -l\ . . . , l' (23) 

"1 J Sphere 

where 



«n2 



,1/ 



R 



A n0 (r) gr 3 dr 



~ J*{A n2 {r) + 3B n2 (r)) or 3 dr 



(24.a) 



(24.b) 



The functions A n i(r), B n [{r) are given in Appendix B, and R is the sphere's radius. To our reassur- 
ance, only the monopole and quadrupole sphere modes survive, as seen by the presence of the factors 
Sin m fl23|)- The final series is thus a relatively simple one, even in spite of its generality^]: 



oo oo 
u ( x ^) = 2^ — ""OOWffnO W + 1^ ~ 



n=l Un0 
where, it is recalled, 



n=l 



J2 u n2m (x)^ 2 2 m) (t) 



(t > 0) (25) 



9nl 



{lm \t) 



g^ lm \t') S muj nl (t-t')dt' , (1 = 0,2; 



m 



(26) 



Equation (|25| ) constitutes the sphere's response to an arbitrary tidal GW perturbation, and will be 
used to analyse the sensitivity of the spherical detector in the next section. Before doing so, however, 
a few comments on the antenna's signal deconvolution capabilities, within the context of a completely 
general metric theory of GWs, are in order. 



3.1 The deconvolution problem 

Let us first of all take the Fourier transform of (25): 

/oo 
u(x, t) e~ 
-oo 

This is seen to be 



■iujt 



(It 



(27) 



CXJ 

U(x,w) = -J2—"nOo(x)G(° \u;)[5(u;-u Jn o)-5(u J +u; n o)] + 



n=l 



oo 

vr x - a n 2 



l\ y-v U n2 



n=l 



]T u n2m (x)G( 2m )(^) 



m=-2 



(28) 



3 From now on I will drop the label P, meaning spheroidal mode, to ease the notation since toroidal modes no longer 
appear in the formulae. 
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where G^ lm \uj) are the Fourier transforms of g( lm \t), respectively: 

G^ lm \u)= g {lm) {t) e~ iujt dt (29) 



o 



The 5-function factors are of course idealisations corresponding to infinitely long integration times 
and infinitely narrow resonance line-widths — but the essentials of the ensuing discussion will not be 
affected by those idealisations. 

If the measuring system were (ideally) sensitive to all frequencies, filters could be applied to examine 
the antenna's oscillations at each monopole and quadrupole frequency: a single transducer would suffice 
to reveal G^ 00 \uj) around the monopole frequencies uj n o, whilst five (placed at suitable positions) would 
be required to calculate the five degenerate amplitudes G^ 2m \oj) around the quadrupole frequencies 
uj n 2- Once the six functions G^ lm \uj) would have thus been determined, inverse Fourier transforms 
would give us the functions g^ lm \t), and thereby the six Riemann tensor components Roioj(t) through 



inversion of the second of equations (|20|), i.e., as an expansion like (15) — only with g's instead of 5's. 
Deconvolution would then be complete. 

Well, not quite. . . Knowledge of the Riemann tensor in the laboratory frame coordinates is not 
really sufficient to say the waveform has been completely deconvolved, unless we also know the source 
position in the sky. There clearly are two possibilities: 

i) The source position is known ahead of time by some other astronomical observation methods. 
Let me rush to emphasise that, far from trivial or uninteresting, this is a very important case to 
consider, specially during the first stages of GW Astronomy, when any reported GW event will 
have to be thoroughly checked by all possible means. 

If the incidence direction is known, then a rotation must be applied to the just obtained quantities 
RoiOj(t), which takes the laboratory z-axis into coincidence with the incoming wave propagation 
vector. A classification procedure must thereafter be applied to the so transformed Riemann 
tensor in order to see which is the theory (or class of theories) compatible with the actual ob- 



servations. Such classification procedure has been described in detail in [20]; see also [49 for an 
updated discussion. 

The spherical antenna is thus seen to have the capability of furnishing the analyst sufficient 
information to discern amongst different competing theories of GW physics, whenever the wave 
incidence direction is known prior to detection. 

ii) The source position is not known at detection time. This makes things more complex, since the 
above rotation between the laboratory and GW frames cannot be performed. 

In order to deconvolve the incidence direction in this case, a specific theory of the GWs must be 
assumed — a given choice being made on the basis of whatever prior information is available or, 



simply, dictated by the the decision to probe a particular theory. Wagoner and Paik |6l[ propose 
a method which is useful both for GR and BD theory, their idea being simple and elegant at the 
same time: since neither of these theories predicts the excitation of the m=±l quadrupole modes 
of the wave, the source position is determined precisely by the rotation angles which, when applied 
to the laboratory axes, cause the amplitudes of those antenna modes to vanish; the rotated frame 
is thereby associated to the GW natural frame. 

A generalisation of this idea can conceivably be found on the basis of a detailed — and possibly 
rather casuistic — analysis of the canonical forms of of the Riemann tensor for a list of theories 
of gravity, along the following line of argument: any one particular theory will be characterised 
by certain (homogeneous) canonical relationships amongst the monopole and quadrupole com- 
ponents of the Riemann tensor, g^ lm ^(t), and so enforcement of those relations upon rotation of 
the laboratory frame axes should enable determination of the rotation angles or, equivalently, of 
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the incoming radiation incidence direction. Scalar-tensor theories e.g. have g^ 2:izl \t) = in their 
canonical forms, hence Wagoner and Paik's proposal for this particular case. 

Before any deconvolution procedure is triggered, however, it is very important to make sure that 
it will be viable. More precisely, since the transformation from the laboratory to the ultimate 
canonical frame is going to be linear, invariants must be preserved. This means that, even if the 
source position is unknown, certain theories will forthrightly be vetoed by the observed Roioj(t) if 
their predicted invariants are incompatible with the observed ones. To give but an easy example, 
if Roioj(t) is observed to have a non-null trace RoiQi(t), then a veto on GR will be readily served, 
and therefore no algorithm based on that theory should be applied [pi]] . 

I would like to make a final remark here. Assume a direction deconvolution procedure has been 
successfully carried through to the end on the basis of certain GW theory, so that the analyst 
comes up with a pair of numbers (0, ip) expressing the source's coordinates in the sky. Of course, 
these numbers will represent the actual source position only if the assumed theory is correct. 
Now, how do we know it is correct? Strictly speaking, "correctness" of a scientific theory is 
an asymptotic concept — in the sense that the possibility always remains open that new facts 
be eventually discovered which contradict the theory — , and so reliability of the estimate (9, ip) 
of the source position can only be assessed in practice in terms of the consistency between the 
assumed theory and whatever experimental evidence is available to date, including, indeed, GW 
measurements themselves. It is thus very important to have a method to verify that the estimate 
(9, <p) does not contradict the theory which enabled its very determination. 

Such verification is a logical absurdity if only one measurement of position is available; this 
happens for instance if the recorded signal is a short burst of radiation, and so two antennas 
are at least necessary to check consistency in that case. The test would proceed as a check that 
the time delay between reception of the signal at both detectors is consistent with the calculated 
(9, ip*f\, given their relative position and the wave propagation speed predicted by the assumed 
theory. If, on the other hand, the signal being tracked is a long duration signal, then a single 
antenna may be sufficient to perform the test by looking at the observed Doppler patterns and 
checking them against those expected with the given (9,ip). 

The above considerations have been made ignoring noise in the detector and monitor systems. A 
fundamental constraint introduced by noise is that it makes the antenna bandwidth limited in sensitivity. 
As a consequence, any deconvolution procedure is deemed to be incomplete or, rather, ambiguous ||, 
since information about the signal can possibly be retrieved only within a reduced bandwidth, whilst 
the rest will be lost. I thus come to a detailed discussion of the sensitivity of the spherical GW antenna 
in the next section. 

4 The sensitivity parameters 

I will consider successively amplitude and energy sensitivities; the first leads to the concept of transfer 
function, while the second to that of absorption cross section. I devote separate subsections to analyse 
each of them in some detail. 

4.1 The transfer function 

A widely used and useful concept in linear system theory is that of transfer function [^] . It is defined 
as the Fourier transform of the system's impulse response, or as the system's impedance/admittance, 
and can be inferred from the frequency response function (f28|), 

4 Note that the two detectors will agree on the same (6,ip), even if the assumed theory is wrong, since the sphere 
deformations will be the same if caused by the same signal. 
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We recall from the previous section that the sphere is a multimode device — due to its monopole 
and five-fold degenerate quadrupole modes. It is expedient to define a multimode or vector transfer 
function as a useful construct which encompasses all six different modes into a single conceptual block, 
according to 

U(x, w) = z(Q) ( x > w ) ° {a) M ( 3 °) 

a 

where G^ a \uj) are the six driving terms G^ lm \uj) given in (|29|). The transfer function is Z( q )(x,cj), 
and its "vector" character alluded above is reflected by the multimode index a. Looking at (|28j) it is 
readily seen that 

Z (00) (x,u;) = - V — u n00 (x) [5(uj-uj n0 ) - 5(uj+uj n0 )] (31.a) 

00 

Z( 2m )(x, W ) = - V ^ Un2m (x) [5(uj-uj n2 ) - 5(uj+uj n2 )] (m = -2,...,2) (31.b) 
* fi ^2 

As we observe in these formulae, the sphere's sensitivity to monopole excitations is governed by 
o-no/^no, and to quadrupole ones by a n2 /uj n2 . Closed expressions happen to exist for a n o and a n2 ; using 
the notation of Appendix B, they are 



OnO 3C(w,0) j 2 (q n oR) 

8tt q n0 R 



a n 2 _ 3C(n, 2) 
i? 8vr 



(33(k n2 R) 3 - — f3 1 {q n2 R) — — — 

Qn2tt &n2 K n2 K 



(32.a) 



(32.b) 



Numerical investigation of the behaviour of these coefficients shows that they decay asymptotically 

as n~ 2 : 

a n i — > const x n {66) 

Likewise, it is found that the frequencies uj h q and uj n2 diverge like n for large n, so that Z( Q )(x,w) 
drops as w -3 for large uj. Figures 6 and 7 display a symbolic plot of u> 3 Z(°°)(x, uj) and w 3 Z( 2m )(x,tj), 
respectively, which illustrates the situation: monopole modes soon reach the asymptotic regime, while 
there appear to be 3 sub-families of quadrupole modes regularly intertwined; the asymptotic regime 
for these sub-families is more irregularly reached. Note also the perfectly regular alternate changes of 
phase (by tt radians) in both monopole and each quadrupole family. 

The sharp fall in sensitivity of a sphere for higher frequency modes (n -3 ) indicates that only the 
lowest ones stand a chance of being observable in an actual GW antenna. I report in Table I the 
numerical values of the relevant parameters for the first few monopole and quadrupole modes. The 
reason for the last (fourth) columns will become clear later. 

4.2 The absorption cross section 

Let us calculate now the energy of the oscillating sphere. We first define the spectral energy density at 
frequency uj, which is naturally given byQ 

J T is the integration time — assumed very large. The peaks in the 5-functions diverge like T/w. 
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Table 1: First few monopole (left) and quadrupole (right) sphere parameters, for a cr=0.33 material. 
First and second columns on either side of the central line number the modes and give the corresponding 
eigenvalue; rows are intertwined in order of ascending frequency, which is proportional to kR — see 
( p^) below. Third columns contain the a n Q and a n i coefficients defined in equations (24); the fourth 
columns display the cross section ratios (fciodio/^noOno) 2 and (^12^12 / 'k^a^) 2 for higher frequency 
modes, respectively, taking as reference the lowest in each family — cf. equations (^ 



n 



k n oR 



1 5.432 0.214 



2 12.138 -3.772xl(T 2 6.46 



3 18.492 1.600xl(T 2 15.49 



n 


k n2 R 


0-112/ R 


o~nlo- n i 


1 


2.650 


0.328 


1 


2 


5.088 


0.106 


2.61 


3 


8.617 


-1.907 xl0~ 2 


27.95 


4 


10.917 


-9.101xl0~ 3 


76.42 


5 


12.280 


1.387xl0~ 2 


25.99 


6 


15.347 


6.879xl0~ 3 


67.87 



W(u) = ^ I \u 2 |U(x,^)| 2 gd 3 x 
1 J Solid 2 



and can be easily evaluated: 

-. 00 1 2 

W(u) = -ttMY j E E a li Gilm) ^) 



(34) 



(35) 



n=l 1=0,2 m=-l 



The energy at any one spectral frequency uj n i is obtained by integration of the spectral density in 
a narrow interval around uj = ±oj n f 



E(uj n i 



27T 



In this case, 



E{u nl )= l -Ma 2 nl E |G (/m) C 



Un2, 



m=—l 



1 = 0,2 



(36) 



(37) 



The sensitivity parameter associated with the vibrational energy of the modes is the detector's 
absorption cross section, defined as the energy it absorbs per unit incident GW spectral flux density, 
or 



0"abs {V) 



E(uj) 
${uj) 



(38) 



where <£(cj) is the number of joules per square metre and Hz carried by the GW at frequency uj as it 
passes by the antenna. Thus, for the frequencies of interest, 
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1 G(°°)Ko) 2 

CTabs(Wno) = ^a n0 — — r — (39. a) 



1 2 E™=- 2 G( 2 ™)(u, n2 ) 2 , 

CTabs(Wn2) = 2 Ma «2 $£7^) ^ ' 

These quantities have very precise values, but such values can only be calculated on the basis 
of a specific underlying theory of the GW physics. In the absence of such theory, neither <3?(u;) nor 
G" m >{u) can possibly be calculated, since they are not theory independent quantities. To date, only 
GR calculations have been reported in the literature [13, p3l |6l[]. As I will now show, even though the 



fractions in the rhs of (|39| ) are not theory independent, some very general results can still be obtained 
about the sphere's cross section within the context of metric theories of the gravitational interaction. 
To do so, it will be necessary to go into a short digression on the general nature of weak metric GWs. 

No matter which is the (metric) theory which happens to be the "correct one" to describe gravitation, 
it is beyond reasonable doubt that any GWs reaching the Earth ought to be very weak. The linear 
approximation should therefore be an extremely good one to describe the propagating field variables 
in the neighbourhood of the detector. In such circumstances, the field equations can be derived from a 
Poincare invariant variational principle based on an action integral of the type 

C(iJA,^A,n)d 4 x (40) 

where the Lagrangian density £ is a quadratic functional of the field variables iPa(x) and their 
space-time derivatives tpA,fi(x); these variables include the metric perturbations h^ v , plus any other 
fields required by the specific theory under consideration — e.g. a scalar field in the theory of Brans- 
Dicke, etc. The requirement that C be quadratic ensures that the Euler-Lagrange equations of motion 
are linear. 

The energy and momentum transported by the waves can be calculated in this formalism in terms 
of the components t^ v of the canonical energy-momentum tensor^] 



(x '^?^- £ ' r (41) 

The flux energy density, or Poynting, vector is given by Si = c 2 r° l , i.e., 

BC 

S(x,t)=c 3 £-^W>4 (42) 

where '= d/dt. Any GW hitting the antenna will be seen plane, due to the enormous distance to the 
source. If k is the incidence direction (normal to the wave front), then the fields will depend on the 
variable ct — k x, so that the GW energy reaching the detector per unit time and area is 

BC 

<P(t) = k-S(x,t) = -c 2 £ ttT ^ (43) 

A d VA 

where x is the sphere's centre position relative to the source — which is fixed, and so its dependence 



can be safely dropped in the lhs of the above expression. The important thing to note in equation (43) 



6 This tensor is not symmetric in general, but can be symmetrized by a standard method due to Belinfante B3|. For 
the considerations which follow in this paper it is unnecessary to go into those details, and the canonical form Ell) will 
be sufficient. 
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is that it tells us that cj>{t) can be written as a quadratic form in the time derivatives of the fields 
As a consequence, the spectral density 5>(o>), defined by 



j_ oo 4>(t)dt = j o *(«)^ (44) 



can be ascertained to factorise as 



= uj 2 $ (lj) (45) 



where < l ) o( w ) is again a quadratic function of the Fourier transforms *$>a(^) of the fields ipA- On the 
other hand, the functions G^ lm \uj) in (|39| ) which, it is recalled, are the Fourier transforms of g( lm \t) 
in (|20|), contain second order derivatives of the metric fields /i^, and therefore of a/Z the fields ipA as a 
result of the theory's field equations. Since we are considering plane wave solutions to those equations, 



all derivatives can be reduced to time derivatives — just like in (43) above. We can thus write 



G {lm \uj) = -uj 2 ^ lm \u) (46) 

with fy( lm \uj) suitable linear combinations of the ^a(^)- Replacing the last two equations into ([39]) 
and manipulating dimensions expediently, we come to the remarkable result that 



Cabs l^nOj 



GMv\ 



3 (knO^no) 2 



GMv? 



(k n 2a n 2Y 



(47.a) 
(47.b) 



where v 2 



(2+2cr) _1 v 2 , v s being the speed of sound in the detector's material, and a its Poisson ratio; 
G is the Gravitational constant. The "remarkable" about the above is that the coefficients Ks(&) 
and Kq(R) are independent of frequency: they exclusively depend on the underlying gravitation theory, 
which I symbolically denote by To see that this is the case, it is enough to consider a monochromatic 
incident wave: since the coefficients Ks{^) and JTq(N) happen to be invariant with respect to field 
amplitude settlings, this means they will only depend on the amplitudes' relative weights, i.e., on the 
field equations' specific structure. 

By way of example, it is interesting to see what the results for General Relativity (GR) and Brans- 
Dicke (BD) theory are. After somehow lengthy algebra it is found that 








16 7T 2 

15 



(48) 



and 



i + 



N = BD 



(3+2H) 2 



(49) 



16 7T 2 
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1 + i (3 + 20)- 2 k 



1 + 



hU 



In the latter formulae, f2 is the usual Brans-Dicke parameter lo [11 1, renamed here to avoid confusion 
with frequency, and A: is a dimensionless parameter, generally of order one, depending on the source's 
properties [36]. As is well known, GR is obtained in the limit f2 — > oo of BD fl6"2|| ; the quoted results 
are of course in agreement with that limit. 

Incidentally, an interesting consequence of the above equations is that the presence of a scalar field 
in the theory of Brans and Dicke causes not only the monopole sphere's modes to be excited, but also the 
m=0 quadrupole ones; what we see in equations ( |49] ) is that precisely 5/6 of the total energy extracted 
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from the scalar wave goes into the antenna's monopole modes, whilst there is still a remaining 1/6 
which is communicated to the quadrupoles, independently of the values of Q and 

This somehow non-intuitive result finds its explanation in the structure of the Riemann tensor in 
BD theory, in which the excess Roioj with respect to General Relativity happens not to be proportional 

to the scalar part E^ ^ , but to a combination of E^ °^ and E^ ^ . More in-depth analysis of these facts 
has been investigated in reference [[10(1 . 

Equations (47) show that, no matter which is the gravity theory assumed, the sphere's absorption 
cross sections for higher modes scale as the successive coefficients (k n oa n o) 2 and (/c n 2«n2) 2 for monopole 
and quadrupole modes, respectively. In particular, the result quoted in [13] that cross section for the 
second quadrupole mode is 2.61 times less than that for the first, assuming GR, is in fact valid, as we 
now see, independently of which is the (metric) theory of gravity actually governing GW physics. The 
fourth columns in Table 4.1 display these scaling properties. It is seen that the drop in cross section 
from the first to the second monopole mode is as high as 6.46. It should however be stressed that 
the frequency of such mode would be over 4 kHz for a (likely) sphere whose fundamental quadrupole 
frequency be 900 Hz [13]. Note finally the asymptotic cross section drop as n~ 2 for large n — cf. 
equation (33) and the ensuing paragraph. 



5 Conclusion 

The main purpose of this paper has been to set up a sound mathematical formalism to address with as 
much generality as possible any questions related to the interaction between a resonant antenna and a 
weak incoming GW, with much special emphasis on the homogeneous sphere. New results have been 
found along this line, such as the scaling properties of cross sections for higher frequency modes, or the 
sensitivity of the antenna to arbitrary metric GWs; also, new ideas have been put forward regarding the 
direction deconvolution problem within the context of an arbitrary metric theory of GW physics. Less 
spectacularly, the full machinery has also been applied to produce independent checks of previously 
published results. 

The whole investigation reported herein has been developed with no a priori assumptions about 
any specific (metric) theory of the GWs, and is therefore very general. "Too general solutions" are 
often impractical in science; here, however, the "very general" appears to be rather "cheap", as seen 
in the results expressed by the equations of section 3 above. An immediate consequence is that solid 
elastic detectors of GWs (and, in particular, spheres) offer, as a matter of principle, the possibility 
of probing any given theory of GW physics with just as much effort as it would take, e.g., to probe 
General Relativity: the vector transfer function of section 4 supplies the requisite theoretical vehicle 
for the purpose. 

An important question, however, has not been considered in this paper. This is the transducer 
problem: the sphere's oscillations can only be revealed to the observer by means of suitable (usually 
electromechanical) transducers. These devices, however, are not neutral, i.e., they couple to the an- 
tenna's motions, thereby exercising a back action on it which must be taken into consideration if one 
is to correctly interpret the system's readout. Preliminary studies and proposals have already been 
published ]^0| , but further work is clearly needed for a more thorough understanding of the problems 
involved. 

Progress in this direction is currently being made — which I expect to report on shortly. The 
formalism developed in this paper provides basic support to that further workP[ 

7 Note however that since monopole and quadrupole detector modes occur at different frequencies, this particular 
distribution of energy may not be seen if the sphere's vibrations are monitored at a single resonance. 

8 This was underway in 1995, and is now complete. The results are presented right below, from page ^ on. 
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Appendices 

A.l Algebra of the ^-matrices 

Let e x ,ey,e z be three orthonormal Cartesian vectors defining the sphere's laboratory reference frame. 
We define the equivalent triad 

e (0) =e 2 , = ^(e x ±ie y ) (50) 

having the properties 

e *K) . e M = 5m , m t m, m ' = -1,0,1 (51) 

We say that the vectors ( |50|) are the natural basis for the 1=1 irreducible representation of the 
rotation group; they behave under arbitrary rotations precisely like the spherical harmonics Yi m (n). In 
particular, if a rotation of angle a around the z-axis is applied to the original frame then 

e (±D ^exp(±ia)e( ±1 ) , -» (52) 

Higher rank tensors have specific multipole characteristics depending on the number of tensor in- 
dices, and the above basis lends itself to reveal those characteristics, too. For example, the five dimen- 
sional linear space of traceless symmetric tensors supports the 1=2 irreducible representation of the 
rotation group, while a tensor's trace is an invariant. A general symmetric tensor can be expressed as 
an "orthogonal" sum of a traceless symmetric tensor and a multiple of the unit tensor. A convenient 
basis to expand any such tensor is the following: 



e (o) e (i) + e« g> e<°) 

e (i) ® e (-!) + e^ 1 ) $ 



e (o) g, e (-i) + e (-i) 

-2e(°) ®e(°) 
e (i) + e (°) ® e(°) 



,(0) 



(53.a) 
(53.b) 
(53.c) 
(53.d) 



The elements ( |53.a ) get multiplied by e^ 2ta in a rotation of angle a around the z-axis, respectively, 
the ( pTb|) by e ±ia , and (px|) and fl53j| ) are invariant, as is readily seen. These properties define the 
"spin characteristics" of the corresponding tensors. Also, the five elements (|53.a|) ~( [53.cj ) are traceless 
tensors, while ( 53. d| ) is the unit tensor. Any symmetric tensor can be expressed as a linear combination 
of the six (53), and the respective coefficients carry the information about the weights of the different 
monopole and quadrupole components of the tensor. 

Equations (|0]) in the text are the matrix representation of the above tensors in the Cartesian basis 
e^e^e^ except that they are multiplied by suitable coefficients to ensure that the conditions 



E 



(lm) 



Y lm (n) , J = 0,2; 



m 



(54) 



'% 3 n i n 3 

where n = x/|x| is the radial unit vector, hold. They are arbitrary, but expedient for the calculations 
in this paper. The following orthogonality relations can be easily established: 
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p *(2m') „(2m) 15 r F *(00) p (2m) _ n F *(00) p (00) _ 3 /rrn 

with the indices m,m! running from —2 to 2, and with an understood sum over the repeated % and j. 
It is also easy to prove the closure properties 

4 ( ° 0) < ] + | E 4 (2m) ^ = I W fl + 6»6 jk ) (56) 

m=— 2 



Equations (|55|) and ( p6|) constitute the completeness equations of the S-matrix basis of Euclidean 
symmetric tensors. 

A. 2 The sphere's spectrum and wave-functions 

This Appendix is intended to give a rather complete summary of the frequency spectrum and eigenmodes 



of a uniform elastic sphere. Although this is a classical problem in Elasticity Theory [42], some of the 



results which follow have never been published so far. Also, its scope is to serve as reference for notation, 



etc., in future work — see ensuing article in page 25. 

The uniform^ elastic sphere's normal modes are obtained as the solutions to the eigenvalue equation 

/iV 2 u+ (A + m) V(V-u) = -u?q\i (57) 

with the boundary conditions that its surface be free of any tensions and/or tractions; this is expressed 
by the equations [32] 

Oij rij = at r = R (58) 
where R is the sphere's radius, n the outward normal, and o~ij the stress tensor 

Oij = A u kk 5ij + 2[i mj (59) 

with Uij = \[uij + Uj^i), the strain tensor, and the Lame coefficients |3^]. 

Like any differentiable vector field, u(x) can be expressed as a sum of an irrotational vector and a 
divergence-free vector, 

u(x) = U irrot .(x) + U div _ free (x) , (60) 

say; on substituting this into equation (|57|), and after a few easy manipulations, one can see that 

(V 2 + k 2 ) u div _ frce (x) = , (V 2 + q 2 ) u irrot . (x) = (61) 

where 

/x A + 2/Lt 

Now the irrotational component can generically be expressed as the gradient of a scalar function, 

i.e., 

Uirrot. (x) = V0(x) (63) 

while there are two linearly independent divergence-free components which, as can be readily verified, 
are 



9 By uniform I mean its density g is constant throughout the solid in the unperturbed state. 
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u 



(1) 

div— free 



x)=L^(x), and ug|_ free (x) = VxL^(x) 



(64) 



where L = -ixxV is the "angular momentum" operator, cf. |21| ], and and ip^ are also scalar 
functions. If ( |63| ) and (|64j) are now respectively substituted in (|6lD, it is found that </>(x), ^(^(x), and 
■0( 2 )(x) satisfy Helmholtz equations: 



(V 2 + A; 2 )V(x) =0 



(V 2 + (7 2 )</)(x)=0 



where ip(x.) stands for either ^^(x) or ip^(x.). Therefore 



<K X ) = ji(qr)Y lm (n) 



^(x) = j;(fcr)y Zm (n) 



(65) 



(66) 



in order to ensure regularity at the centre of the sphere, r=0. Here, ^ is a spherical Bessel function 
— see HI for general conventions on these functions — , and Yi m a spherical harmonic [pTf. Finally thus, 



u(x) = % V<«x) + ^LV>(x) + § V / L. ixl 



where Co,Ci,C 2 are three constants which will be determined by the boundary conditions ( |58|) (the 
denominators under them have been included for notational convenience). After lengthy algebra, those 
conditions can be expressed as the following system of linear equations: 



where 



02{qR) ~ ^ q 2 R 2 (3 (qR)\ C - 1(1 + 1) 0x{kR) C 2 

(3i (kR) C\ 



fafaR) Co - [§ (3 2 (kR) + {^±±1 - l} o (kR) 



There are clearly two families of solutions to fl68| ) : 
i) Toroidal modes. These are characterised by 



Co 



5. a 



(68.b) 
(68.c) 

(69) 



p 1 (kR)=0, C = C 2 = 



(70) 



The frequencies of these modes are independent of A, and thence independent of the material's 
Poisson ratio. Their amplitudes are 



u nim(x) = T nl (r) iLY lm (n) 



(71) 



with 



T n i{r) = C 1 (n,l)ji(k n ir) 



(72) 



and C\(n,l) a dimensionless normalisation constant determined by the general formula (JsJ) ; k n iR 
is the n-th root of the first equation (170) for a given I. 



18 



ii) Spheroidal modes. These correspond to 



det 



(73) 



f3 2 (qR)-±q 2 R 2 f3 (qR) 1(1 + 1) ^(kR) 

/3 1 (qR) \ 2 (kR) + - l} Po(kR) 

and C\ = 0. The frequencies of these modes do depend on the Poisson ratio, and their amplitudes 



arc 



u nZm( x ) = A m(r) Yi m (n) n - B n i(r) inx LYi m (n) 
where A n i(r) and B n i(r) have the somewhat complicated form 



(74) 



A n i{r) 



C(n,l) 



B n i(r) = C(n,l) 



HknlR)j'l(qnir) - 1(1 + 1) ^ Pl(qnlR) Jj ^^ 
r (u n\ MVnir) q nl {k nl r ji(k nl r)}' 

P3{KnlR) T Pl{qnlR> 7 

qnl r Knl Knl r 



with accents denoting derivatives with respect to implied (dimensionless) arguments, 



p 3 ( z ) = lf3 2 (z) + { l -^-l} /3 (z) 



(75.a) 
(75.b) 



(76) 



and C(n,l) a new normalisation constant. It is understood that q n \ and k n \ are obtained after 
the (transcendental) equation ( f73[ ) has been solved for to — cf. equation (|62|). 



In actual practice equations ( |70| ) and ( |73| ) are solved for the dimensionless quantity kR, which 
will hereafter be called the eigenvalue. In view of (|62|), the relationship between the latter and the 
measurable frequencies (in Hz) is given by 



v = 



uJ 

2^ 



V2 kR 
qR?) 2vr 



(77) 



It is more useful to express the frequencies in terms of the Poisson ratio, a, and of the speed of 
sound v s in the selected material. For this the following formulas are required — see e.g. [j32f: 



where Y is the Young modulus, related to the Lame coefficients and the Poisson ratio by 

(3A + 2/^)// _ A 



Y 



A + /i 



Hence, 



2(1 + 0-)//, a 



(kR) v s 



2(A + /i) 



(78) 



(79) 



(80) 



2vr^/TTo : R 

Equation (|80| ) provides a suitable transformation formula from abstract number eigenvalues (kR) 
into physical frequencies v, for given material's properties and sizes. 

Tables ||] and ^ respectively display a set of values of (kR) for toroidal and spheroidal modes. While 
GWs can only couple to quadrupole and monopole modes, it is important to have some detailed knowl- 
edge of analytical results, as the sphere's frequency spectrum is rather involved. It often happens, both 
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Table 2: List of a few spheroidal eigenvalues, ordered in columns of ascending harmonics for each 
multipole value. Spheroidal eigenvalues depend on the sphere's material Poisson ratio — although this 
dependence is weak. In this table, values are given for a = 0.33. Note that the table contains all 
eigenvalues less than or equal to 11.024 yet is not exhaustive for values larger than that one; this would 
require to stretch the table horizontally beyond I = 10 — see Figure I for a qualitative inspection of 
trends in eigenvalue progressions. 



n 


I = 


/ = 1 


I = 2 


I = 3 


I = 4 


I = 5 


I = 6 


I = 7 


I = 8 


I = 9 


/ = 10 


1 


5.4322 


3.5895 


2.6497 


3.9489 


5.0662 


6.1118 


7.1223 


8.1129 


9.0909 


10.061 


11.024 


2 


12.138 


7.2306 


5.0878 


6.6959 


8.2994 


9.8529 


11.340 


12.757 


14.111 


15.410 


16.665 


3 


18.492 


8.4906 


8.6168 


9.9720 


11.324 


12.686 


14.066 


15.462 


16.867 


18.272 


19.664 


4 


24.785 


10.728 


10.917 


12.900 


14.467 


15.879 


17.243 


18.589 


19.930 


21.272 


22.619 


5 


31.055 


13.882 


12.280 


14.073 


16.125 


18.159 


19.997 


21.594 


23.043 


24.426 


25.778 



in numerical simulations and in experimental determinations, that it is very difficult to disentangle the 
wealth of observed frequency lines, and to correctly associate them with the corresponding eigenmode. 
Complications are further enhanced by partial degeneracy lifting found in practice (due to broken sym- 
metries), which result in even more frequency lines in the spectrum. Accurate analytic results should 
therefore be very helpful to assist in frequency identification tasks. 



Table 3: List of a few toroidal eigenvalues, ordered in columns of ascending harmonics for each multipole 
value. Unlike spheroidal eigenvalues, toroidal eigenvalues are independent of the sphere's material 
Poisson ratio. Note that the table contains all eigenvalues less than or equal to 12.866 yet is not 
exhaustive for values larger than that one; this would require to stretch the table horizontally beyond 
/ = 11 — see Figure II for a qualitative inspection of trends in eigenvalue progressions. 



n 


I = 1 


I = 2 


I = 3 


I = 4 


I = 5 


/ = 6 


/ = 7 


/ = 8 


I = 9 


/ = 10 


I = 11 


1 


5.7635 


2.5011 


3.8647 


5.0946 


6.2658 


7.4026 


8.599 


9.6210 


10.711 


11.792 


12.866 


2 


9.0950 


7.1360 


8.4449 


9.7125 


10.951 


12.166 


13.365 


14.548 


15.720 


16.882 


18.035 


3 


12.323 


10.515 


11.882 


13.211 


14.511 


15.788 


17.045 


18.287 


19.515 


20.731 


21.937 


4 


15.515 


13.772 


15.175 


16.544 


17.886 


19.204 


20.503 


21.786 


23.055 


24.310 


25.555 


5 


18.689 


16.983 


18.412 


19.809 


21.181 


22.530 


23.860 


25.174 


26.473 


27.760 


29.035 



In Figures 1 and 2 a symbolic line diagramme of the two families of frequencies of the sphere's 
spectrum is presented. Spheroidal eigenvalues have been plotted for the Poisson ratio <r=0.33. Although 
only the 1=0 and 1=2 spheroidal series couple to GW tidal forces, the plots include other eigenvalues, 
as they can be useful both in bench experiments — cf. equation ( |l2|) above — and for vetoing purposes 
in a spherical antenna. 

Figures 3, 4 and 5 contain plots of the first three monopole and quadrupole functions T n i(r), A n i(r) 
and B n i(r), always for <r=0.33. T n o(r) and B n o(r) have however been omitted; this is because they are 



multiplied by an identically zero angular coefficient in the amplitude formulae (71) and (|74|), Indeed, 
monopole vibrations are spherically symmetric, i.e., purely radial. 
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Figure I The homogeneous sphere spheroidal eigenvalues for a few multipole families. Only the 
1=0 and 1=2 families couple to metric GWs, so the rest are given for completeness and non-directly-GW 
uses. Note that there are fewer monopole than any other /-pole modes. The lowest frequency is the 
first quadrupole. The diagram corresponds to a sphere with Poisson ratio cr=0.33. Frequencies can be 
obtained from the plotted values through equation (|8QD for any specific case. 



Figure II The homogeneous sphere toroidal eigenvalues. None of these couple to GWs, but 
knowledge of them can be useful for vetoing purposes. These eigenvalues are independent of the 
material's Poisson ratio. To obtain actual frequencies from plotted values, use (|S0[). The lowest toroidal 
eigenvalue is kR = 2.5011, with 1=2, and happens to be the absolute minimum sphere's eigenvalue. 
Compared to the spheroidal kR = 2.6497, also with 1=2, its frequency is 5.61% smaller. Note also that 
there are no monopole toroidal modes. 



Figure III First three spheroidal monopole radial functions A n o(r) (n = 1,2,3), equation (75. a). 



Figure IV First three spheroidal quadrupole radial functions A n 2(r) (continuous line) and B n 2(r) 
(broken line) (n = 1,2,3), equations (ffq). 



Figure V First three toroidal quadrupole radial functions T n 2{r) (n = 1,2,3), equation ([721). A 
common feature to these radial functions (also in the two previous Figures) is that they present a nodal 
point at the origin (r = 0), while the sphere's surface (r/R = 1) has a non-zero amplitude value, which 
is largest (in absolute value) for the lowest n in each group. 



Figure VI The scalar component Z^ 00 ^(x, uj) of the multimode transfer function, ( 31.a| ). The 



diagram actually displays uj 3 Z*- 00 ) (x, uj), so asymptotic behaviours are better appreciated. It is given in 
units of n/gR, and a factor (vr/i) u n oo( x )j the eigenmode amplitude, has been omitted, too. <5-function 



amplitudes are symbolically taken as 1. Note that the asymptotic regime, given by equation (33), is 
quickly reached. 



Figure VII The quadrupole component Z^ 2 " 1 ) (x, uj) of the multimode transfer function, ( |31.b ). 



The same prescriptions of Figure 6 apply here; the plot is therefore independent of the value of m. 
Note the presence of three sub-families of peaks; asymptotic regimes are reached with variable speed 
for these sub-families, and less rapidly than for monopole modes, anyway. 
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Abstract 

Apart from omnidirectional, a solid elastic sphere is a natural multi-mode and multi-frequency 
device for the detection of Gravitational Waves (GW). Motion sensing in a spherical GW detector 
thus requires a multiple set of transducers attached to it at suitable locations. If these are resonant 
then they exert a significant back action on the larger sphere and, as a consequence, the joint 
dynamics of the entire system must be properly understood before reliable conclusions can be drawn 
from its readout. In this paper, I present and develop an analytic approach to study such dynamics 
which generalises currently existing ones and clarifies their actual range of validity. In addition, the 
new formalism shows that there actually exist resonator layouts alternative to the highly symmetric 
TIG A, potentially having interesting properties. One of these (I will call it PHC), which only requires 
five resonators per quadrupole mode sensed, and has mode channels, will be described in detail. Also, 
the perturbative nature of the proposed approach makes it very well adapted to systematically assess 
the consequences of realistic mistunings in the device parameters by robust analytic methods. In 
order to check the real value of the mathematical model, its predictions have been confronted with 
experimental data from the LSU prototype detector TIGA, and agreement between both is found 
to consistently reach a satisfactory precision of four decimal places. 



1 Introduction 

The idea of using a solid elastic sphere as a gravitational wave (GW) antenna is almost as old as that 
of using cylindrical bars: as far back as 1971 Forward published a paper [p2j | in which he assessed 
some of the potentialities offered by a spherical solid for that purpose. It was however Weber's ongoing 
philosophy and practice of using bars which eventually prevailed and developed up to the present date, 
with the highly sophisticated and sensitive ultra-cryogenic systems currently in operation — see |l 



and [16] for reviews and bibliography. With few exceptions [Q, |6lj ], spherical detectors fell into oblivion 
for years, but interest in them strongly re-emerged in the early 1990's, and an important number of 
research articles have been published since which address a wide variety of problems in GW spherical 
detector science. At the same time, international collaboration has intensified, and prospects for the 
actual construction of large spherical GW observatories (in the range of ~100 tons) are being currently 
considered in several countries 0, even in a variant hollow shape [|l7]] . 

A spherical antenna is obviously omnidirectional but, most important, it is also a natural multi- 
mode device, i.e., when suitably monitored, it can generate information on all the GW amplitudes and 
incidence direction p3{| , a capability which possesses no other individual GW detector, whether resonant 
or interferometric [jyj. Furthermore, a spherical antenna could also reveal the eventual existence of 



monopole gravitational radiation, or set thresholds on it [10]. 



The theoretical explanation of these facts is to be found in the unique matching between the GW 



amplitude structure and that of the sphere oscillation eigenmodes [39, 41]: a general metric GW gener- 
ates a tidal field of forces in an elastic body which is given in terms of the "electric" components Roioj (t) 
of the Riemann tensor at its centre of mass by the following formula, see equation (]i~9| ) above 

10 There are collaborations in Brazil, Holland, Italy and Spain. 
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f GW (x,t) = Yl f (Zm) (x)5 ( ' m) (t) (81) 

1=0 and 2 
m=— 

where f (x) are "tidal form factors" , while g^ m ^ (t) are specific linear combinations of the Riemann 
tensor components Roioj(t) which carry all the dynamical information on the GW's monopole (1 = 0) 
and quadrupole (I = 2) amplitudes. It is precisely these amplitudes, g( lm ^ (t), which a GW detector aims 
to measure. 

On the other hand, a free elastic sphere has two families of oscillation eigenmodes, so called toroidal 
and spheroidal modes, and modes within either family group into ascending series of /-pole harmonics, 
each of whose frequencies is (2Z+l)-fold degenerate — see |37j for full details. It so happens that 
only monopole and/or quadrupole spheroidal modes can possibly be excited by an incoming metric 
GW |J, and their GW driven amplitudes are directly proportional to the wave amplitudes g( lm \t) 
of equation (|8l|). It is this very fact which makes of the spherical detector such a natural one for 



GW observations [37|. In addition, a spherical antenna has a significantly higher absorption cross 



section than a cylinder of like fundamental frequency, and also presents good sensitivity at the second 



quadrupole harmonic [13]. 



In order to monitor the GW induced deformations of the sphere motion sensors are required. In 
cylindrical bars, current state of the art technology is based upon resonant transducers 0, p4[ . A res- 
onant transducer consists in a small (compared to the bar) mechanical device possessing a resonance 
frequency accurately tuned to that of the cylinder. This frequency matching causes back-and-forth res- 
onant energy transfer between the two bodies (bar and resonator), which results in turn in mechanically 
amplified oscillations of the smaller resonator. The philosophy of using resonators for motion sensing 
is directly transplantable to a spherical detector — only a multiple set rather than a single resonator is 
required if its potential capabilities as a multi-mode system are to be exploited to satisfaction. 

A most immediate question in a multiple motion sensor system is: where should the sensors be? 
The answer to this basic question naturally depends on design and purpose criteria. Merkowitz and 
Johnson (M&J) made a very appealing proposal consisting in a set of 6 identical resonators coupling 
to the radial motions of the sphere's surface, and occupying the positions of the centres of the 6 non- 
parallel pentagonal faces of a truncated icosahedron [30, 46]. One of the most remarkable properties 



of such layout is that there exist 5 linear combinations of the resonators' readouts which are directly 
proportional to the 5 quadrupole GW amplitudes g( 2m "> (t) of equation (|8l|). M&J call these combinations 
mode channels, and they therefore play a fundamental role in GW signal deconvolution in a real, 



noisy system [49, pi] . In addition, a reduced scale prototype antenna — called TIG A, for Truncated 



Icosahedron Gravitational Antenna — was constructed at Louisiana State University, and its working 



experimentally put to test [45]. The remarkable success of this experiment in almost every detail [47, 
|48l , 50] stands as a vivid proof of the practical feasibility of a spherical GW detector || [52]| . 

Despite its success, the theoretical model proposed by M&J to describe the system dynamics is 
based upon a simplifying assumption that the resonators only couple to to the quadrupole vibration 
modes of the sphere ^(| . While this is seen a posteriori of experimental measurements to be a very 



good approximation [45, 48 1, a deeper physical reason which explains why this happens is missing so 
far. The original motivation for the research I present in this article was to develop a more general 
approach, based on first principles, for the analysis of the resonator problem, very much in the spirit 
of the methodology and results of reference |J7|; this, I thought, would not only provide the necessary 
tools for a rigorous analysis of the system dynamics, but also contribute to improve our understanding 
of the physics of the spherical GW detector. 

Pursuing this programme, I succeeded in setting up and solving the equations of motion for the 
coupled system of sphere plus resonators. The most important characteristic of the solution is that it 
is expressible as a perturbative series expansion in ascending powers of the small parameter r/ 1 / 2 , where 
r] is the ratio between the average resonator's mass and the sphere's mass. The dominant (lowest) 
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order terms in this expansion appear to exactly reproduce Merkowitz and Johnson's equations [46], 
whence a quantitative assessment of their degree of accuracy, as well as of the range of validity of their 
underlying hypotheses obtains; if further precision is required then a well defined procedure for going 
to next (higher) order terms is unambiguously prescribed by the system equations. 

Beyond this, though, the simple and elegant algebra which emerges out of the general scheme has 
enabled the exploration of different resonator layouts, alternative to the unique TIG A of M&J. In 
particular, I found one (3^, 4C] requiring 5 rather than 6 resonators per quadrupole mode sensed and 
possessing the remarkable property that mode channels can be constructed from the system readouts, 
i.e., five linear combinations of the latter which are directly proportional to the five quadrupole GW 
amplitudes. I called this distribution PHC — see below for full details. 

The intrinsically perturbative nature of the proposed approach makes it also particularly well 
adapted to assess the consequences of small defects in the system structure, such as for example sym- 
metry breaking due to suspension attachments, small resonator mistunings and mislocations, etc. This 
has been successfully applied to account for the reported frequency measurements of the LSU TIGA 
prototype [45], which was diametrically drilled for suspension purposes; in particular, discrepancies be- 
tween measured and calculated values (generally affecting only the fourth decimal place) are precisely 
of the theoretically predicted order of magnitude. 

The method has also been applied to analyse the stability of the spherical detector to several 
mistuned parameters, with the result that it is not very sensitive to small construction errors. This 
conforms again to experimental reports [50|, but has the advantage that the argument depends on 
analytic mathematical work rather than on computer simulated procedures — see e.g. |5(| or [58]. 

The paper is structured as follows: in section ^, I present the main physical hypotheses of the model, 
and the general equations of motion. In section |3] a Green function approach to solve those equations is 
set up, and in section f| it is used to assess the system response to both monopole and quadrupole GW 



signals. In section 
channels. Section 



5] I describe in detail the PHC layout, including its frequency spectrum and mode 
contains a few brief considerations on the system response to a hammer stroke 
calibration signal, and finally in section [?] I assess how the different parameter mistunings affect the 
detector's behaviour. The paper closes with a summary of conclusions, and three appendices where the 
heavier mathematical details are made precise for the interested reader. 



2 General equations 

With minor improvements, I shall use the notation of references p^] and | |38| |, some of which is now 
briefly recalled. Consider a solid sphere of mass A4, radius R, (uniform) density g, and elastic Lame 
coefficients A and /x, endowed with a set of J resonators of masses M a and resonance frequencies Q a 
(a = l,. . . ,J), respectively. I shall model the latter as point masses attached to one end of a linear 
spring, whose other end is rigidly linked to the sphere at locations x a — see Figure |2[ The system 
degrees of freedom are given by the field of elastic displacements u(x, t) of the sphere plus the discrete 
set of resonator spring deformations z a (t); equations of motion need to be written down for them, of 
course, and this is my next concern in this section. 

I shall assume that the resonators only move radially, and also that Classical Elasticity theory [32] 
is sufficiently accurate for the present purposes^. In these circumstances we thus have 



d 2 



u 



//V 2 u+ (A + /x)V(V-u) +f(x,t) (82.a) 



• dt 2 

z a (t) = -n 2 a [z a (t)-u a (t)]+C xtcrn&l (t) , a=l,...,J (82.b) 



1 We clearly do not expect relativistic motions in extremely small displacements at typical frequencies in the range of 
1 kHz. 
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Figure 1: Schematic diagram of the coupling model between a solid sphere and a resonator. The 
notation is that in the text, but sub- indices have been dropped for clarity. The dashed-dotted arc 
line on the left indicates the position of the undeformed sphere's surface, and the solid arc its actual 
position. 



where n a = x a /i? is the outward pointing normal at the the a-th resonator's attachment point, and 



u a {t) = n a -u(x a ,i) , a = l,...,J (83) 

is the radial deformation of the sphere's surface at x a . A dot (') is an abbreviation for time derivative. 
The term in square brackets in ( |S2.b ) is thus the spring deformation — q(t) in Figure ||. 

f (x, t) in the rhs of ( |32.a| ) contains the density of all non-internal forces acting on the sphere, which 
is expediently split into a component due the resonators' back action and an external action proper, 
which can be a GW signal, a calibration signal, etc. Then 



f (X, t) = frcsonators (x, t) + f ex ternal (x, t ) (84) 

Finally, ^ Brtornal (t) in the rhs of $5ty is the force per unit mass (acceleration) acting on the a-th 
resonator due to external agents. 

Given the hypothesis that the resonators are point masses, the following holds: 

J 

fre S onators(x, t) = [ z a{t) ~ U a (t)\ <5 (3) (x - X a ) n a (85) 

a=l 

where 5^ is the three dimensional Dirac density function. 

The external forces I shall be considering in this paper will be gravitational wave signals, and also 
a simple calibration signal, a perpendicular hammer stroke. GW driving terms, c.f. equation (|sTl), can 
be written 

f GW (x,t) = f( 00 )(x)< 7 ( °)(t) + £ f( 2 ™)(x)«/ 2m )(t) (86) 

m=-2 

for a general metric wave — see |57j for explicit formulas and technical details. While the spatial 
coefficients f (x) are pure form factors associated to the tidal character of a GW excitation, it is the 
time dependent factors g(' m ) (t) which carry the specific information on the incoming GW. The purpose 
of a GW detector is to determine the latter coefficients on the basis of suitable measurements. 
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If a GW sweeps the observatory then the resonators themselves will also be affected, of course. 
They will be driven, relative to the sphere's centre, by a tidal acceleration which, since they only move 
radially, is given by 

£ w (t) = c 2 R 0i0j {t) x a:i n aJ , a = 1, . . . , J (87) 

where RoiOj(t) are the "electric" components of the GW Riemann tensor at the centre of the sphere. 
These can be easily manipulated to give0 

£ W (t) = R £ Y lm (n a )g {lm Ht) , a = l,...,J (88) 

1=0 and 2 

where R is the sphere's radius. 

I shall also be later considering the response of the system to a particular calibration signal, consist- 
ing in a hammer stroke with intensity f"o, delivered perpendicularly to the sphere's surface at point xo: 

fstroke(x, t) = f*o ^ 3 )(x-x )(5(t) (89) 

which is modeled as an impulsive force in both space and time variables. Unlike GW tides, a hammer 
stroke will be applied on the sphere's surface, so it has no direct effect on the resonators. In other 
words, 

e° ke W = 0, a = l,...,J (90) 
The fundamental equations thus finally read: 



d 2 u 

= /^V 2 u + (A + ^)V(V-u) + 

J 

M b n 2 b [z b (t) - u b (t)} 5®(x - x 6 ) n b + f cxtcrnal (x, t) (91.a) 

6=1 

z a (t) = -n 2 a [z a (t)-u a (t)}+C xter ^\t) , a = l,...,J (91.b) 

where f e xternai( x ; t) will be given by either (|86| ) or (|39l), as the case may be. Likewise, £f ternal (£) will be 
given by ( |88| ) or (p0|), respectively. The remainder of this paper will be concerned with finding solutions 
to the system of coupled differential equations (pip, and with their meaning and consequences. 

3 Green function formalism 



In order to solve equations (|9l|) I shall resort to Green function formalism. The essentials of this 
procedure in the context of the present problem can be found in detail in reference 37]; more specific 



technicalities are given in appendix AT| 



By means of such formalism equations (^) become the following integro-differential system: 

u a {t) = < xternal (t) + E% f : K ab (t-t')[z b (t')-u b (t')]dt' (92.a) 
Zait) = C xm \t) - n 2 a [za(t) - u a (t)] , a = 1, . . . , J (92.b) 



12 Yim(n) are spherical harmonics |21| — see also the multipole expansion of Roioj(t) in reference |37 
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where it| xteraal (i) = n a • u external (x a , i) , and u external ( x , t) is the bare (i.e., without attached resonators) 
sphere's response to the external forces f ex ternai( x ) t) in the rhs of (|91.a[ ). K ab (t) is a kernel matrix 
defined by the following weighted sum of diadic products of wavefunctionsP|: 

K ab {t) = ttlY^ ^n 1 [n6-u^(x 6 )] [n a -UAr(x a )] smuj N t (93) 

JV 

Finally, the mass ratios of the resonators to the entire sphere are defined by 

rib = ^ , b = 1, . . . , J (94) 

and will be small parameters in a real device. 

Before proceeding further, let us briefly pause for a qualitative inspection of equations (P3). Equa- 



tion ( 92.a ) shows that the sphere's surface deformations u a (t) are made up of two contributions: one 
due to the action of external agents (GWs or other), contained in it| xternal (i), and another one due to 
coupling to the resonators. The latter is commanded by the small parameters r) b , and correlates to all 
of the sphere's spheroidal eigenmodes through the kernel matrix K ab (t). This has consequences for GW 
detectors, for even though GWs may only couple to quadrupole and monopole^ spheroidal modes of 
the free sphere |?7|, [|, attachment of resonators causes, as we see, coupling between these and the other 
modes of the antenna; conversely, the latter back-act on the former, too. As I shall shortly prove, such 
undesirable effects can be minimised by suitably tuning the resonators' frequencies. 



3.1 Laplace transform domain equations 

A solution to equations (|9^) will now be attempted. Equation ( |92.a| ) is an integral equation belonging in 
the general class of Volterra equations ||5S| ], but the usual iterative solution to it by repeated substitution 
of u b (t) into the kernel integral is not viable here due to the dynamical contribution of z b (t), which is 
in turn governed by the differential equation ( |92.b ). 

A better suited method to solve this integro- differential system is to Laplace-transform it. I denote 
the Laplace transform of a generic function of time f(t) with a caret (") on its symbol, e.g., 

rod 

f(s)= / f(t)e~ st dt (95) 
Jo 

and make the assumption that the system is at rest before an instant of time, t = 0, say, or 

u(x, 0) = u(x, 0) = z a (0) = z a (0) = (96) 
Equations ( |92"|) then adopt the equivalent form 



J 

u a (s) = uT tei ™\s)-Y,VbK ab ( S )[z b (s)-u b ( S )} (97.a) 

6=1 

s 2 B a (s) = £ xternal (s) - n 2 a [z a (s) - u a (s)} , a = 1, . . . , J (97.b) 

for which use has been made of the convolution theorem for Laplace transforms^. A further simplifi- 
cation is accomplished if we consider that we shall in practice be only concerned with the measurable 
quantities 

13 The capitalised index N will often be used to imply the multiple index {nlm} which characterises the sphere's 
wave-functions. 

14 Monopole modes only exist in scalar-tensor theories of gravity, such as e.g. Brans-Dicke Jnj|; General Relativity of 
course does not belong in this category. 

15 This theorem states, it is recalled, that the Laplace transform of the convolution product of two functions is the 
arithmetic product of their respective Laplace transforms. 
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q a {t) = z a (t) - u a (t) , a=l,...,J (98) 

representing the resonators' actual elastic deformations — cf. Figure ^. It is readily seen that these 
verify the following: 



J 

E 

6=1 



S 2 



o~ab + rib 2 Q2 K ab (s) 

S + 5 l n 



„2 ^external I \ 

M = - 7 r^T2 u7 m \s) + ^T^T - « = 1, • • • , J (99) 



Equations (|99|) constitute a significant simplification of the original problem, as they are a set of 
just J algebraic rather than integral or differential equations. We must solve them for the unknowns 
q a (s), then perform inverse Laplace transforms to revert to q a (t). I come to this next. 

4 System response to a Gravitational Wave 

Our concern now is the actual system response when it is acted upon by an incoming GW. It will be 
calculated by making a number of simplifying assumptions, more precisely: 

) The detector is perfectly spherical. 

) The resonators have identical masses and resonance frequencies. 

) The resonators' frequency is accurately matched to one of the sphere's oscillation eigenfrequencies. 

It will be shown below (section [?]) that a real system can be appropriately treated as one which 
deviates by definite amounts from this idealised construct. Therefore detailed knowledge of the ideal 
system behaviour is essential for all purposes: such is the justification for the above simplifications. 



The wave- functions u n ; m (x) of an elastic sphere can be found in reference [37] in full detail, and I 



shall keep the notation of that paper for them. The Laplace transform of the kernel matrix (|93|) can 



thus be expressed as — see equation ( |158|) in appendix [A. 1 



nl S +U) nl 47r n l S + ^nl 

where the last term simply defines the quantities Xat>^ ■ Note that the sums here extend over the entire 
spectrum of the solid sphere. 

The assumption that all the resonators are identical simply means that 

m = ... = Vj = rj , n 1= ...=nj = n (101) 

The third hypothesis makes reference to the fundamental idea behind using resonators, which is to 
have them tuned to one of the frequencies of the sphere's spectrum. This is expressed by 

n = uj nolo (102) 

where u) no i is a specific and fixed frequency of the spheroidal spectrum. 

In a GW detector it will only make sense to choose Iq = or Iq = 2, as only monopole and quadrupole 
sphere modes couple to the incoming signal; in practice, no will refer to the first, or perhaps second 
harmonic |l3j . I shall however keep the generic expression (|102| ) for the time being in order to encompass 
all the possibilities with a unified notation. 

Based on the above hypotheses, equation (^) can be rewritten in the form 
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J 

E 

6=1 



2„2 



(nZ) 

2 \ ^afe 



where i s ) ls the Laplace transform of 



i.e., 



-GW 



+ n 2 



(s) + 



e w (*) 

s 2 + ft 2 



(O 



4 GW (*) 



i? 2 r/ m K)5 (/m) (s) 



w„ / ) 

(103) 
(104) 



2=0 and 2 
m=— i,...,Z 



As mentioned at the end of the previous section, the matrix in the lhs of ( |103| ) must now be inverted; 
this will give us an expression for q a (s), whose inverse Laplace transform will take us back to the time 
domain. A simple glance at the equation suffices however to grasp the unsurmountable difficulties of 
accomplishing this analytically. 

Thankfully, though, a perturbative approach is applicable when the masses of the resonators are 
small compared to the mass of the whole sphere, i.e., when the inequality 



7] < 1 (105) 

holds. I shall henceforth assume that this is the case, as also is with cylindrical bar resonant transducers. 
It is shown in appendix [A.2] that the perturbative series happens in ascending powers of n 1 ^ 2 , rather 
than n itself, and that the lowest order contribution has the form 

?«(*) = V~ 1/2 £ A a lm) (s; n) (a) + O(0) , a = 1, . . . , J (106) 

l,m 

where O(0) stands for terms of order rp or smaller. Here, Aa (s; O) is a transfer function matrix which 
relates linearly the system response q a ( s ) to the GW amplitudes g^ lm \s), in the usual sense that q a (t) 
is given by the convolution product of the signal g^ m ) (t) with the time domain expression, 
of A. a lm \s;Q). The detector is thus seen to act as a linear filter on the GW signal, whose frequency 
response is characterised by the properties of A a lrn \s; $7). More specifically, the filter has a number of 
characteristic frequencies which correspond to the imaginary parts of the poles of Ai im '(s; Q). As also 
shown in appendix |A.2| , these frequencies are the symmetric pairs 



-tf ( 1±\Ft^ \A no i (R)\ CaV 1/2 I +0( V ) , a = l,...,J (107) 



47T 

where £ 2 is the a-th eigenvalue of the Legendre matrix 



P lo (n a -n b ) , a, 6 = 1,..., J (108) 

associated to the multipole (Zo) selected for tuning — see ( |102[) . These frequency pairs correspond to 
beats, typical of resonantly coupled oscillating systems — we shall find them again in section | in a 
particularly illuminating example. 

Equation ( |106| ) neatly displays the amplification coefficient r/ -1 / 2 of the resonators' motion ampli- 
tudes, which corresponds to the familiar resonant energy transfer in coupled systems of linear oscilla- 
tors £§. 

The specific form of the transfer function matrix Aq (s; Q) depends both on the selected mode to 
tune the resonator frequency O and on the resonator distribution geometry. I now come to a discussion 
of these. 
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4.1 Monopole gravitational radiation sensing 

General Relativity, as is well known, forbids monopole GW radiation. More general metric theories, 
e.g. Brans-Dicke do however predict this kind of radiation. It appears that a spherical antenna 
is potentially sensitive to monopole waves, so it can serve the purpose of thresholding, or eventually 
detecting them. It is therefore relevant to consider the system response to scalar waves. 

This clearly requires that the resonator set be tuned to a monopole harmonic of the sphere, i.e., 



Q = ca 



nO 



(*o = 0) 



(109) 



where n tags the chosen harmonic — most likely the first (n = 1) in a thinkable device. 
Since Pq(z) = 1 (for all z) the eigenvalues of Po( n a - nf,) are, clearly, 



Ci 2 = J , cl 



c5 = o 



(110) 



for any resonator distribution. The tuned mode frequency thus splits into a single strongly coupled 
pair: 



J 



c4 = ^ | 1 ± \l — \A n0 (R)\ rj 1 ^ ) + 0(7,) , n = uj n0 



(111) 



The A-matrix of equation (106) is seen to be in this case 



Ai /m) ( S ;^„o) 
whence the system response is 



QnO 1 

77 2 



(s 2 + <4) _1 - (s 2 + JL) 



'm0 



:n2 



-1/2 ("1) J 1 



2 , 2 



2 , 2 
S + U!_ 



9 m (s)+O(0) , 



1,...,J (113) 



regardless of resonator positions. The overlap coefficient a n Q is given by ( 24. a ), and can be calculated 
by means of numerical computer programmes. By way of example, a\o/R = 0.214, and 020/ R = —0.038 
for the first two harmonics. 

A few interesting facts are displayed by equation (|113| ). First, as already stressed, it is seen that if 
the resonators are tuned to a monopole detector frequency then only monopole wave amplitudes couple 
strongly to the system, even if quadrupole radiation amplitudes are significantly high at the observation 
frequencies uj±. Also, the amplitudes q a (s) are equal for all a, as corresponds to the spherical symmetry 
of monopole sphere's oscillations, and are proportional to J -1 / 2 , a factor we should indeed expect as 
an indication that GW energy is evenly distributed amongst all the resonators. A single transducer 
suffices to experimentally determine the only monopole GW amplitude g^ 00 \s), of course, but (|113| ) 
provides the system response if more than one sensor is mounted on the antenna for whatever reasons. 



4.2 Quadrupole gravitational radiation sensing 

I now consider the more interesting case of quadrupole motion sensing. The choice is now, clearly, 



n = cu n2 , (Z = 2) (114) 

where n labels the chosen harmonic — most likely the first (n = 1) or the second (n = 2) in a practical 
system. The evaluation of the A-matrix is now considerably more involved |37]], yet a remarkably 
elegant form is found for it: 
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A(M( S ;u, n2 ) = (-1)" 7f a n2 E (e | f(? + - (* 2 + - 2 ^ 



6=1 U=^° 



(c)„,(c)* 



>2m( n b)^2 

(115) 

where v a is the c-th normalised eigenvector of i-2(n a -nft), associated to the non-null eigenvalue £ 2 . Let 



me stress that equation (115) explicitly shows that at most 5 pairs of modes, of frequencies w c ±, couple 
strongly to quadrupole GW amplitudes, no matter how many resonators in excess of 5 are mounted on 
the sphere. The tidal overlap coefficients a2 n can also be calculated using (p4b|) , and give for the first 
two harmonics 



Q12 
R 



0.328 , 



The system response is thus 



Q22 
R 



0.106 



(116) 



la(s) 



-1/2 




2 i 2 
S +UJ r4 



-1 



2 i 2 



-i 



(c) (c)* 



E l2m(n 6 )5 (2m) (s) + O(0) , a = l...,J 



(117) 



Equation ( |117| ) is completely general, i.e., it is valid for any resonator configuration over the sphere's 
surface, and for any number of resonators. It describes precisely how all 5 GW amplitudes g^ 2m ^(s) 
interact with all 5 strongly coupled system modes; like before, only quadrupole wave amplitudes are seen 
in the detector (to leading order) when 0, = oj n 2 , even if the incoming wave carries significant monopole 
energy at the frequencies u> c ±. 

The degree of generality and algebraic simplicity of ( |117| ) is new in the literature. As we shall now 
see, it makes possible a systematic search for different resonator distributions and their properties. 



5 The PHC configuration 

Merkowitz and Johnson's TIGA [^] is highly symmetric, and is the minimal set with maximum de- 
generacy, i.e., all the non-null eigenvalues C, a are equal. To accomplish this, however, 6 rather than 5 
resonators are required on the sphere's surface. Since there are just 5 quadrupole GW amplitudes one 



may wonder whether there are alternative layouts with only 5 resonators. Equation (117) is completely 



general, so it can be searched for an answer to this question. In reference [38] a specific proposal was 
made, which I now describe in detail. 

In pursuing a search for 5 resonator sets I found that distributions having a sphere diameter as an 
axis of pentagonal symmetrical exhibit a rather appealing structure. More specifically, let the resonators 
be located at the spherical positions 

2tt 

6 a = a (alia), ip a = (a - 1) — , a = l,...,5 (118) 

5 

The eigenvalues and eigenvectors of P2(n a -rife) are easily calculated: 

Co 2 = | (3 cos 2 a -if , Ci = C-i = f sin 2 a cos 2 a , C 2 2 = C 2 2 = f sm 4 a (119.a) 
v { r ] = Jf C- 1 Y 2m (n a ) , m = -2,...,2, o = l,...,5 (119.b) 



By this I mean resonators are placed along a parallel of the sphere every 72° 
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so the A-matrix is also considerably simple in structure in this case: 



A^(s;u; n2 ) = 
where the notation 



f4jF - L l 

^ a n2 W 2 



Y 2m (n a )5 l2 , PHC (120) 



2 

UJ m± 



1 ± 



m 



-2,. ..,2 



fl21 



has been used. As seen in these formulas, the /i«e expected pairs of frequencies actually reduce to three, 
so pentagonal distributions keep a certain degree of degeneracy, too. The most important distinguishing 
characteristic of the general pentagonal layout is best displayed by the explicit system response: 



Qa(s) 



a n 2 
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+ 



{^[(^o 2 +r-(* 2 + 
'(* 2 +-i+r-( s2+ ^-) 
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-1 



^2oK)5 (20) (s) 

Y 22 (n a ) gW(s) + y 2 - 2 (n a ) g^ 2 ~ 2 \s) 



(122) 



This equation indicates that different wave amplitudes selectively couple to different detector fre- 
quencies. This should be considered a very remarkable fact, for it thence follows that simple inspection 
of the system readout spectrum^ immediately reveals whether a given wave amplitude g 2m (s) is present 
in the incoming signal or not. 

Pentagonal configurations also admit mode channels, which are easily constructed from ( |122| ) thanks 
to the orthonormality property of the eigenvectors ( |119.b| ) : 



y {m) (s) 



5 

E 

0=1 



(m)i 



Qa(s) 



f] 1/2 a n2 - 



-1 



g {2m) (s) + O(0) 



(123) 



These are almost identical to the TIG A mode channels Ji6|], the only difference being that each 
mode channel comes now at a single specific frequency pair uj m ±. 

Mode channels are fundamental in signal deconvolution algorithms in noisy systems |4^, pl[ . Pentag- 
onal resonator configurations should thus be considered non-trivial candidates for a real GW detector. 

Based on these facts one may next ask which is a suitable transducer distribution with an axis of 
pentagonal symmetry. Figure [2] shows a plot of the eigenvalues ( 119. a| ) as functions of a, the angular 
distance of the resonator set from the symmetry axis. Several criteria may be adopted to select a 
specific choice in view of this graph. An interesting one can be arrived at by the following argument. If 
for ease of mounting, stability, etc., it is desirable to have the detector milled into a close-to-spherical 
polyhedric shape0 then polyhedra with axes of pentagonal symmetry must be searched. The number 
of quasi regular convex polyhedra is of course finite — there actually are only 18 of them [^8[ — , 
and I found a particularly appealing one in the so called pentagonal hexacontahedron (PHC), displayed 
in Figure ||, left — see also |J8|. This is a 60 face polyhedron, whose faces are the identical irregular 
pentagons of Figure 0, right. The PHC admits an inscribed sphere which is tangent to each face at the 



17 In a noiseless system, of course. 

18 This is the philosophy suggested and experimentally implemented by Merkowitz and Johnson at LSU. 
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Figure 2: The three distinct eigenvalues £ m (m = 0,1,2) as functions of the distance of the resonator 
parallel's co-latitude a relative to the axis of symmetry of the distribution, cf. equation (|119.a|) . 

central point marked in the Figure. It is clearly to this point that a resonator should be attached so as 
to simulate an as perfect as possible spherical distribution. 

The PHC is considerably spherical: the ratio of its volume to that of the inscribed sphere is 1.057, 
which quite favourably compares to the value of 1.153 for the ratio of the circumscribed sphere to the 
TI volume. If the frequency pairs u> m ± are now requested to be as evenly spaced as possible, compatible 
with the PHC face orientations, then the choice a = 67.617° is unambiguously singled out. Hence 



uj 0± = u 12 (l ±0.5756 r? 1/2 ) , wi± = un (l ± 0.8787 ?? 1/2 ) , u 2 ± = W12 (l ± 1-0668 -q 



1/2 



(124) 



for instance for f2 = u>i2, the first quadrupole harmonic. Figure || shows this frequency spectrum together 
with the multiply degenerate TIGA for comparison. 

The criterion leading to the PHC proposal is of course not unique, and alternatives can be consid- 
ered. For example, if the 5 faces of a regular icosahedron are selected for sensor mounting (a = 63.45°) 
then a four-fold degenerate pair plus a single non-degenerate pair is obtained; if the resonator parallel 
is 50° or 22.6° away from the "north pole" then the three frequencies wo+, and uj 2 + are equally 
spaced; etc. The number of choices is virtually infinite if the sphere is not milled into a polyhedric 
shape |H ||]. 



Let me finally recall that the complete PHC proposal [38] was made with the idea of building an 
as complete as possible spherical GW antenna, which amounts to making it sensitive at the first two 
quadrupole frequencies and at the first monopole one. This would take advantage of the good sphere 
cross section at the second quadrupole harmonic |13[| , and would enable measuring (or thresholding) 
eventual monopole GW radiation. Now, the system pattern matrix Aa m \s;Q) has identical structure 
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Figure 3: To the left, the pentagonal hexacontahedron shape. Certain faces are marked to indicate 
resonator positions in a specific proposal — see text — as follows: a square for resonators tuned to the 
first quadrupole frequency, a triangle for the second, and a star for the monopole. On the right we see 
the (pentagonal) face of the polyhedron. A few details about it: the confluence point of the dotted lines 
at the centre is the tangency point of the inscribed sphere to the PHC; the labeled angles have values 
a = 61.863°, [3 = 87.205°; the angles at the T-vertices are all equal, and their value is 118.1366°, while 
the angle at P is 67.4536°; the ratio of a long edge (e.g. PT\) to a short one (e.g. T1T2) is 1.74985, 
and the radius of the inscribed sphere is twice the long edge of the pentagon, R = 2PT\. 



37 



-1.2 -0.8 -0.4 0.0 0.4 0.8 1.2 



TIGA 



o 

03 



PHC 



Frequency 



Figure 4: Compared line spectrum of a coupled TIGA and a PHC resonator layout in an ideally 
spherical system. The weakly coupled central frequency in the TIGA is drawn dashed. The frequency 
pair is 5-fold degenerate for this layout, while the two outer pairs of the PHC are doubly degenerate 
each, and the inner pair is non-degenerate. Units in abscissas are r] l / 2 Vl, and the central value, labeled 
0.0, corresponds to Q. 
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for all the harmonics of a given I series — see ( 112 ) and flllq ) — , and so too identical criteria for resonator 
layout design apply to either set of transducers, respectively tuned to u±2 and 1022- The PHC proposal is 
best described graphically in Figure ^, left: a second set of resonators, tuned to the second quadrupole 
harmonic W22 can be placed in an equivalent position in the "southern hemisphere", and an eleventh 
resonator tuned to the first monopole frequency wio is added at an arbitrary position. It is not difficult 
to see, by the general methods outlined earlier on in this paper, that cross interaction between these 
three sets of resonators is only second order in n 1 ^ 2 , therefore weak. 

A spherical GW detector with such a set of altogether 11 transducers would be a very complete multi- 
mode multi-frequency device with an unprecedented capacity as an individual antenna. Amongst other, 
it would practically enable monitoring of coalescing binary chirp signals by means of a rather robust 
double passage method [15], a prospect which was considered so far possible only with broadband long 
baseline laser interferometers ITq , pi]] , and is almost unthinkable with currently operating cylindrical 
bars. 



6 A calibration signal: hammer stroke 

This section is a brief digression from the main streamline of the paper. I propose to assess now the 
system response to a particular, but useful, calibration signal: a perpendicular hammer stroke. 

Let us first go back to equation ( p9| ) and replace n^ xternal (s) in its rhs with that corresponding to a 
hammer stroke, which is easily calculated — cf. appendix [A.l| : 

-stroke^ = _ £ \ Anl (R)\ 2 PKria-no) , a = 1, . . . , J (125) 

nl S +UJ nl 

where no are the spherical coordinates of the hit point on the sphere, and /o = n o - fo/-^- Clearly, the 
hammer stroke excites all of the sphere's vibration eigenmodes, as it has a completly flat spectrum. 

The coupled system resonances are again those calculated in appendix [A.2| . The same procedures 
described in section || for a GW excitation can now be pursued to obtain 
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6=1 



vi c) vj c) * 



Pi(n b -n ) + O(0) (126) 



where a = 1,..., J, when the system is tuned to the n/-th spheroidal harmonic, i.e., Vt = u) n i. It is 
immediately seen from here that the system response to this signal when the resonators are tuned to a 
monopole frequency is given by 



-1/2 



(-1) 



^= \A n0 (R)\ I 



(s 2 + u,*)' 1 - (s 2 + cu 2 ) 







<^n0 



(127) 



an expression which holds for all a, and is independent of either the resonator layout or the hit point, 
which in particular prevents any determination of the latter, as obviously expected. The frequencies 
u>± are those of ( |111| ), and we find here again a global factor J -1 / 2 , as also expected. 

Consider next the situation when quadrupole tuning is implemented, Q = cj n 2- Only the PHC and 
TIGA configurations will be addressed, as more general cases are not quite as interesting at this point. 
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Figure 5: Simulated response of a PHC to a hammer stroke: the time series and their respective spectra, 
both for direct resonator readouts and mode channels. Note that while the former are not simple beats, 
the latter are. 



40 



6.1 PHC and TIG A response to a hammer stroke 



Expanding equation ( 126 ) by substitution of the eigenvalues Cm and eigenvectors v a 
readily finds that the system response is given by 



(m) 



of the PHC, one 



q a (s)=V- 1,2 k\j^\An2{R)\ £ ~ 

m=— 2 



2 | 2 



with a = 1, ... ,5, and the mode channels by 
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C- 1 Y 2m (n a )Yi m (n ) (128) 



» (m) w 



-V2/0 |^n2(i?)k 



*2m( n o) 



m 



(129) 



These equations indicate that the system response q a (t) is a superposition of three different beatsf^\, 
while the mode channels are single beats each, but with differing modulation frequencies. This is 
represented graphically in Figure [|, where we see the result of a numerical simulation of the PHC 
response to a hammer stroke, delivered to the solid at a given location. The readouts q a (t) are somewhat 
complex time series, whose frequency spectrum shows three pairs of peaks — in fact, the lines in the 
ideal spectrum of Figure ^. The mode channels on the other hand are pure beats, whose spectra consist 
of the individually separate pairs of the just mentioned peaks. 

The response of the TIGA layout to a hammer stroke has been described in detail by Merkowitz 
and Johnson — see e.g. reference ^8|. The present formalism does of course lead to the results obtained 
by them; in the notation of this paper, we have 
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Y 2 * m (n ) , m = -2,...,2 (130.b) 



for the system response and the mode channels, respectively, where 

T 



^2(^)1^/2 +0(7?) , 



1,...,6 



(131) 



are the five-fold degenerate frequency pairs corresponding to the TIGA distribution. Comparison of 
the mode channels shows that they are identical for PHC and TIGA, except that the former come at 
different frequencies depending on the index m. One might perhaps say that the PHC gives rise to 
a sort of "Zeeman splitting" of the TIGA degenerate frequencies, which can be attributed to an axial 
symmetry breaking of that resonator distribution: the PHC mode channels partly split up the otherwise 
degenerate multiplet into its components. 



7 Symmetry defects 

So far we have made the assumption that the sphere is perfectly symmetric, that the resonators are 
identical, that their locations on the sphere's surface are ideally accurate, etc. This is of course un- 
realistic. So I propose to address now how departures from such ideal conditions affect the system 



A beat is a modulated oscillation of the form sin — u)-)t cosQt, where u> + and w_ are nearby frequencies, and 



U+ +LU- —2Q. The Laplace transform of such function of time is precisely (fi/2) 
higher order terms in the difference lo+ — ui- , which in this case is proportional to r\ 



up to 



1/2 
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behaviour. As we shall see, the system is rather robust, in a sense to be made precise shortly, against 
a number of small defects. 

In order to quantitatively assess ideality failures I shall adopt a philosophy which is naturally 
suggested by the results already obtained in an ideal system. It is as follows. 

As seen in previous sections, the solution to the general equations (|9^) must be given as a perturbative 
series expansion in ascending powers of the small quantity r/ 1 / 2 . This is clearly a fact not related to the 
system's symmetries, so it will survive symmetry breakings. It is therefore appropriate to parametrise 
deviations from ideality in terms of suitable powers of ?? 1//2 , in order to address them consistently with 
the order of accuracy of the series solution to the equations of motion. An example will better illustrate 
the situation. 

In a perfectly ideal spherical detector the system frequencies are given by equations ( |107| ). Now, 
if a small departure from e.g. spherical symmetry is present in the system then we expect that a 
correspondingly small correction to those equations will be required. Which specific correction to the 
formula will actually happen can be qualitatively assessed by a consistency argument: if symmetry 
defects are of order n 1 / 2 then equations (107) will be significantly altered in their n 1 / 2 terms; if on 



the other hand such defects are of order rj or smaller then any modifications to equations (107) will 
be swallowed into the O(0) terms, and the more important n 1 / 2 terms will remain unaffected by the 
symmetry failure. One can say in this case that the system is robust against that symmetry breaking. 

More generally, this argument can be extended to see that the only system defects standing a chance 
to have any influences on lowest order ideal system behaviour are defects of order rj 1 ^ 2 relative to an 
ideal configuration. Defects of such order are however not necessarily guaranteed to be significant, and 
a specific analysis is required for each concrete parameter in order to see whether or not the system 
response is robust against the considered parameter deviations. Let us now go into the quantitative 
detail. 

Let P be one of the system parameters, e.g. a sphere frequency, or a resonator mass or location, 
etc. Let -Pideal be the numerical value this parameter has in an ideal detector, and let P rea i be its value 
in the real case. These two will be assumed to differ by terms of order 

t ? 1 /2 ) i.e., 

^Preal = J Pidoal(l+P?? 1/2 ) (132) 



For a given system, p is readily determined adopting ( 132[) as the definition of P r eab once a suitable 



hypothesis has been made as to which is the value of -Pideal- In order for the following procedure to 
make sensible sense it is clearly required that p be of order 1 or, at least, appreciably larger than r/ 1 / 2 . 
Should p thus calculated from (132) happen to be too small, i.e., of order 77 1 / 2 itself or smaller, then 



the system will be considered robust as regards the affected parameter. 
7.1 The suspended sphere 

An earth based observatory obviously requires a suspension mechanism for the large sphere. If a nodal 
point suspension is e.g. selected then a diametral bore has to be drilled across the sphere liq] . The most 
immediate consequence of this is that spherical symmetry is broken, what in turn results in degeneracy 
lifting of the free spectral frequencies uj n i, which now split up into multiplcts u^nim 

(m = -l,...,l). The 

resonators' frequency f2 cannot therefore be matched to the frequency w no ; , but at most to one of 
the 

w no/o m ' Si ^ n this subsection I keep the hypothesis — to be relaxed later, see below — that all the 
resonators are identical, and assume that f2 falls within the span of the multiplet of the w no ; om 's. Then 

<z om = ^ 2 (l+lW /2 ) , m = -l ,...,l (133) 
The coupled frequencies, i.e., the roots of equation (|166|), will now be searched. The kernel matrix 



K a b(s) is however no longer given by (100), due the removed degeneracy of uj n i, and we must stick to 
its general expression (|155|) , or 



42 



Kab(s) = 
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(nlm) 



Following the steps of appendix [A.l| we now seek the roots of the equation 



(134) 
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(135) 

Since relates to uJ no i om through equati on ( | 1 33|) we see that the roots of (135) fall again into either 
of the two categories ( |169| ) (see Appendix [A.2| ), i.e., roots close to ±ifl and roots close to ±iw n / m 
(nl^nolo). I shall exclusively concentrate on the former now. Direct substitution of the series ( |169.a ) 
into ( t!35|) yields the following equation for the coefficient xi : 



det 



i lo 
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(136) 



to which it reduces when p ri 



This is a variation of ([L7 1 

The solutions to (|136|) no longer come in symmetric pairs, like (107) 



0, i.e., when there is full degeneracy. 

Rather, there are 21q+1+J 



of them, with a maximum number of 2(2Zo+l) non-identically zero roots if J> 2/q+iP|. For example, 
if we choose to select the resonators' frequency close to a quadrupole multiplet (lo = 2) then fll36| ) has 
at most 5+J non-null roots, with a maximum ten no matter how many resonators in excess of 5 we 
attach to the sphere. Modes associated to null roots of (|136| ) can be seen to be weakly coupled, just like 
in a free sphere, i.e., their amplitudes are smaller than those of the strongly coupled ones by factors of 
order rf-l 2 . 

In order to assess the reliability of this method I have applied it to see what are its predictions for 
a real system. To this end, data taken with the TIGA prototype at LSI/0 were used to confront with. 
The TIGA was drilled and suspended from the centre, so its first quadrupole frequency split up into a 
multiplet of five frequencies. Their reportedly measured values are 



cjiso = 3249 Hz , W121 = 3238 Hz , lo 12 _i = 3236 Hz , u 122 = 3224 Hz , u 12 -2 = 3223 Hz , (137) 

All 6 resonators were equal, and had the following characteristic frequency and mass, respectively: 

1 



J! = 3241 Hz , rj 
Substituting these values into ( |133[ ) it is seen that 



1762.45 



Po = 0.2075 , p! = -0.0777 , p_i = -0.1036 , p 2 = -0.4393 , p_ 2 



-0.4650 



(138) 



(139) 



Equation ( 136 ) can now be readily solved, once the resonator positions are fed into the matrices 



Xab ■ Such positions correspond to the pentagonal faces of a truncated icosahedron. Merkowitz [45] 



gives a complete account of all the measured system frequencies as resonators are progressively attached 
to the selected faces, beginning with one and ending with six. Figure |6] graphically displays the exper- 
imentally reported frequencies along with those calculated theoretically by solving equation ( |136D . In 
Table 1 I give the numerical values. As can be seen, coincidence between theoretical predictions and 

20 This is a mathematical fact, whose proof is relatively cumbersome, and will be omitted here; let me just mention that 
it has its origin in the linear dependence of more than 2io + l spherical harmonics of order lo. 



These data are contained in reference 



I want to thank Stephen Merkowitz for kindly handing them to me. 
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Table 4: Numerical values of measured and theoretically predicted frequencies (in Hz) for the TIGA 
prototype with varying number of resonators. Relative errors are also shown as parts in 10 4 . The 
calculated values of the tuning and free multiplet frequencies are taken by definition equal to the 
measured ones, and quoted in brackets. In square brackets the frequency of the weakly coupled sixth 
mode in the full, 6 resonator TIGA layout. These data are plotted in Figure |(| 
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Difference 
(parts in 10 4 ) 
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Tuning 


3241 


(3241) 


(0) 


4 reson. 
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3155 


-12 


No reson. 


3223 
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5 reson. 
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experimental data is remarkable: the worst error is 0.2%, while for the most part it is below 0.1%. This 
is a few parts in 10 4 , which is precisely the magnitude of 77, as specified in equation ( 138| ). 



Therefore discrepancies between theoretical predictions and experimental data are exactly as expected, 
i.e., of order 77. In addition, it is also reported in reference ]^8| that the 11-th, weakly coupled mode of 
the TIGA (enclosed in square brackets in Table 1) has a practically zero amplitude, again in excellent 
agreement with the general theoretical predictions about modes beyond the tenth — see paragraph after 
equation (|136| ). 

This is a remarkable result which encouraged a better fit by estimates of next order corrections, i.e., 
Xi of ( 169.a| ). As it turned out, however, matching between theory and experiment does not consistently 



improve in the next step. This is not really that surprising, though, as M&J explicitly state [48] that 
control of the general experimental conditions in which data were obtained had a certain degree of 
tolerance, and they actually show satisfaction that ~1% coincidence between theory and measurement 



is comfortably accomplished. But 1% is two orders of magnitude larger than rj — cf. equation (138) — , 
so failure to refine our frequency estimates to order rj is again fully consistent with the accuracy of 
available real data. 

A word on a technical issue is in order. Merkowitz and Johnson's equations for the TIGA |3(], 46] 
are identical to the equations in this paper to lowest order in rj. Remarkably, though, their reported 
theoretical estimates of the system frequencies are not quite as accurate as those in Table 1 p8| . The 
reason is probably this: in M&J's model these frequencies appear within an algebraic system of 5+J 
linear equations with as many unknowns which has to be solved; here instead the algebraic system has 
only J equations and unknowns, actually equations (|9^). This is a very appreciable difference for the 
range of values of J under consideration. While the roots for the frequencies mathematically coincide in 
both approaches, in actual practice they are estimated, generally by means of computer programmes. 
It is here that problems most likely arise, for the numerical reliability of an algorithm to solve matrix 
equations normally decreases as the rank of the matrix increases. 

7.2 Other mismatched parameters 

We now assess the system sensitivity to small mismatches in resonators' masses, locations and frequen- 
cies. 

7.2.1 Resonator mass mismatches 

If the masses are slightly non-equal then one can write 

M a = V M(l + fi aV l/2 ) , a = l,...,J (140) 
where 77 can be defined e.g. as the ratio of the average resonator mass to the sphere's mass. It is 



immediately obvious from equation (140) that mass non-uniformities of the resonators only affect the 
main equations in second order, since resonator mass non-uniformities result, as we see, in corrections 
of order 77 1 / 2 to 77 1 / 2 itself, which is the very parameter of the perturbative expansions. The system is 



thus clearly robust to mismatches in the resonator masses of the type ( |14 
7.2.2 Errors in resonator locations 

The same happens if the locations of the resonators have tolerances relative to a preselected distribution. 
For let n a be a set of resonator locations, for example the TIGA or the PHC positions, and let n' a be 
the real ones, close to the former: 

to! a = n a + v a n 1 ' 2 , a = l,...,J (141) 
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Figure 6: The frequency spectrum of the TIG A distribution as resonators are progressively added from 
none to 6. Continuous lines correspond to measured values, and dashed lines correspond to their r/ 1 / 2 
theoretical estimates with equation ( |136| ), 
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The values n a determine the eigenvalues £ a in equation ( |107| ), and they also appear as arguments 
to the spherical harmonics in the system response functions of sections |]-^]. It follows from ( |141|) by 
continuity arguments that 



Yi m (n' a ) 

c 



Y lm (n a ) + OirjV 2 ) 



(142.a) 
(142.b) 



Inspection of the equations of sections Q-y shows that both £ a and l/ m (n a ) always appear within 
lowest order terms, and hence that corrections to them of the type (|142j ) will affect those terms in 
second order again. We thus conclude that the system is also robust to small misalignments of the 
resonators relative to pre-established positions. 



7.2.3 Resonator frequency mistunings 

The resonator frequencies may also differ amongst them, so let 

l/2> 



n a = n(i + PaV 1 ' z ) 



I, 



(143) 



To assess the consequences of this, however, we must go back to equation (28) and see what the 

coefficients in its series solutions of the type ( 169.a ) are. The procedure is very similar to that of 

section 7.1, and will not be repeated here; the lowest order coefficient vi is seen to satisfy the algebraic 

1 — 1 2 
equation 



det 



1 



>ab 



— E 

X \ c=0 



xt olo) 5 cb 
X|-Pc 







(144) 



which reduces to ( |170| ) when all the p J s vanish, as expected. This appears to potentially have significant 
effects on our results to lowest order in rj 1 ^ 2 , but a more careful consideration of the facts shows that 
it is probably unrealistic to think of such large tolerances in resonator manufacturing as implied by 
equation ( 143 ) in the first place. In the TIG A experiment, for example [45], an error of order n 1 ^ 2 would 
amount to around 50 Hz of mistuning between resonators, an absurd figure by all means. In a full scale 
sphere (~40 tons, ~3 metres in diameter, ~800 Hz fundamental quadrupole frequency, n ~ 1(T 5 ) the 
same error would amount to between 5 Hz and 10 Hz in resonator mistunings for the lowest frequency. 
This is probably excessive for a capacitive transducer, but may be realistic for an inductive one. With 
this exception, it is thus more appropriate to consider that resonator mistunings are at least of order 
77. If this is the case, though, we see once more that the system is quite insensitive to such mistunings. 

Summing up the results of this section, one can say that the resonator system dynamics is quite 
robust to small (of order 77 1 / 2 ) changes in its various parameters. The important exception is of course 
the effect of suspension drilling, which do result in significant changes relative to the ideally perfect 
device, but which can be relatively easily calculated. The theoretical picture is fully supported by 



experiment, as robustness in the parameters here considered has been reported in the real device [48|. 



8 Conclusions 

A spherical GW antenna is a natural multi-mode device with very rich potential capabilities to detect 
GWs on earth. But such detector is not just a bare sphere, it requires a set of motion sensors to be 
practically useful. It appears that transducers of the resonant type are the best suited ones for an 
efficient performance of the detector. Resonators however significantly interact with the sphere, and 
they affect in particular its frequency spectrum and vibration modes in a specific fashion, which must 
be properly understood before reliable conclusions can be drawn from the system readout. 
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The main objective of this paper has been the construction and development of an elaborate theoret- 
ical model to describe the joint dynamics of a solid elastic sphere and a set of radial motion resonators 
attached to its surface at arbitrary locations, with the purpose to make predictions of the system 
characteristics and response, in principle with arbitrary mathematical precision. 

The solutions to the equations of motion have been shown to be expressible as an ascending series 
in powers of the small "coupling constant" rj, the ratio of the average resonator mass to the mass of the 
larger sphere. The lowest order approximation corresponds to terms of order r/ 1//2 and, to this order, 



previous results [}48| , pq , 44] are recovered. This, I hope, should contribute to clarify the nature of the 
approximations inherent in earlier approaches, and to better understand the physical reason for their 
remarkable accuracy [fHf] . 

In addition, the methods of this paper have permitted us to discover that there can be in fact trans- 
ducer layouts alternative to the highly symmetric TIGA, and having potentially interesting practical 
properties. An example is the PHC distribution, which is based on a pentagonally symmetric set of 5 
rather than 6 resonators per quadruopole mode sensed. This transducer distribution has the property 
that mode channels can be constructed from the resonators' readouts, much in the same way as in the 
TIGA [46|. In the PHC however a new and distinctive characteristic is present: different wave ampli- 
tudes selectively couple to different detector modes having different frequencies, so that the antenna's 
mode channels come at different rather than equal frequencies. The PHC philosophy can be extended 
to make a multi-frequency system by using resonators tuned to the first two quadrupole harmonics of 
the sphere and to the first monopole, an altogether 11 transducer set. 

The assessment of symmetry failure effects, as well as other parameter departures form ideality, 
has also been subjected to analysis. The general scheme is again seen to be very well suited for the 
purpose, as the theory transparently shows that the system is robust against relative disturbances of 
order r/ or smaller in any system parameters, also providing a systematic procedure to assess larger 
tolerances — up to order rj 1 ^ 2 . The system is shown to still be robust to tolerances of this order in 
some of its parameters, whilst it is not to others. Included in the latter group is the effect of spherical 
symmetry breaking due to system suspension in the laboratory, which causes degeneracy lifting of the 
sphere's eigenfrequencies, which split up into multiplets. A strong point is that, by use of mostly 
analytic algorithms, it has been possible to accurately reproduce the reportedly measured frequencies 
of the LSU prototype antenna [^] with the predicted precision of four decimal places. The also 
reported robustness of the system to resonator mislocations ]^8| is too in satisfactory agreement with 
the theoretical predictions. 

The perturbative approach here adopted is naturally open to refined analysis of the system re- 
sponse in higher orders in rj. For example, one can systematically address the weaker coupling of 
non-quadrupole modes, etc. It appears however that such refinements will be largely masked by noise 
in a real system, as shown by Merkowitz and Johnson [ 50 1 , and this must therefore be considered first. 
So the next step is to include noise in the model and see its effect. Stevenson [58] has already made some 
progress in this direction, and partly assessed the characteristics of TIGA and PHC, but more needs 
to be done since not too high signal-to-noise ratios should realistically be be considered in an actual 
GW detector. In particular, mode channels are at the basis of noise correlations and dependencies, as 
well as the errors in GW parameter estimation [51]. I do expect the analytic tools developed in this 
article to provide a powerful framework to address the fundamental problems of noise in a spherical 
GW antenna which, to my knowledge, have not yet received the detailed attention they require. 
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Appendices 

A.l Green functions for the multiple resonator system 

The density of forces in he rhs of equation ( |82.a ) happens to be of the separable type 

f(x,t) = £fM(x) «,(«)(*) (145) 



where a is a suitable label. It is recalled from reference [37] that, in such circumstances, a formal 
solution can be written down for equation ( |32.a| ) in terms of a Green function integral, whereby the 
following orthogonal series expansion obtains: 



u(x, t) = £ £ oft 1 # u*(x) g%> (t) 



(146) 



where 



A*) 
In 



a N 



1 



u^(x) -f( a )(x)d 3 x 



g {a) {t') smu N {t-t')dt' 



(147.a) 



(147.b) 



Here, lon and ujy(x) are the eigenfrequencies and associated normalised wave- functions of the free 
sphere. Also, iV is an abbreviation for a multiple index {nlm}. The generic index a is a label for the 
different pieces of interaction happening in the system. I quote the result of the calculations of the 
terms needed in this paper: 



,(o) 

J resonators, N 



Ma 2 r * I \ 

— il a [n a -u N {x a ) 



Mm') 
JGW,N 



rn 



a = 1, . . . , J 

a n i $w $ mm ' , N = {nlm} , l' = 0, 2 

/stroke,7V = M ~ 1 fo • U^r (x ) 

where the coefficients a n i in ( [148. b ) are overlapping integrals of the type (147. a), and 



I' ...J' 



(148.a) 

(148.b) 
(148.c) 



a (a) (t) 

resonators , iV V" / 



[z a {t') - u a (t)] sin lun (t - t') alt' 



a = 1, . . . , J 



9stroke,N(t) = SmUtft 

If this is replaced into ( |32.a| ) one readily finds 



g {lm) {t') smuj N (t-t')dt' 



u ( x , *) = £ <V u iv(x) | £ -tj fi& [n 6 -u^(x 6 )] ff r ( 2 onatorS)iV (*) + £ /< 

N U=l a 



(149.a) 

(149.b) 
(149.c) 



I external,7V ^external, TV (*) (150) 
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where the label "external" explicitly refers to agents acting upon the system from outside. Two kinds 
of such external actions are considered in this article: those due to GWs and those due to a calibration 
hammer stroke signal. Specifying x = x a in the lhs of (150) and multiply on either side by n a , the 
following is readily found: 



uait) = < xternal W + E% / K ab (t-t')[z b (t')-u b (t')]dt' 

z a (t) = e tcrnal (*) " tf a [ z a (t) - u a (t)] , a = 1, . . . , J 
where rj b = M b /M , < xternal (i) = n a ■ u cxtcrnal (x a , t) , 



u external (x,t) = EE^ 1 / ( 

a N 



i f( a ) i \ (°0 ii\ 

external^ U M X J 5 e xternal,AT \ l ) 



and 



K ab {t) Efl^w^ 1 [n b -u^(x 6 )] [n a -UAr(x Q )] smto N t 



N 



(151.a) 
(151.b) 

(152) 
(153) 



The following bare sphere responses to GWs and hammer strokes (equations! (pl|) and (89), respec- 
tively) can be calculated by direct substitution. The results as best presented as Laplace transform 
domain functions: 



E (E an i Anl f ] ) Y lm (n a ) 9 {lm) (s) - a = l 7 (L>l.a) 



2=0 and 2 \n=l 



uf oke (s) 



E 

III 



fo 



, T \A nl (R)\ z P^K-no) 



a = 1, . . . , J 



(154.b) 



where Yi m ar e spherical harmonics and P[ Legendre polynomials [21|. The calculation of the Laplace 
transform of the kernel matrix (11531) is likewise immediate: 



K ab (s) = E 



«2 



N 



S 2 +UJ% 



[n b -u* N (x b )] [n a -UAr(x a )] 



Given that (see 37] for full details) 

Unim(x) = A ni (r) Yi m (6, <p)n- B n i(r) inx LYj m (0, <p) 



(155) 



(156) 



and that the spheroidal frequencies uj n i are 2/+l-fold degenerate, ( 155 ) can be easily summed over the 
degeneracy index m, to obtain 



K ab {s) = E n » 

nl ° ^^nl 



., i 2~ \ A nl{R)Y 



E Y l*m( n b) Y lm( n a) 
m=—l 



or, equivalently, 



nl S +U nl 
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(157) 



(158) 



where use has been made of the summation formula for the spherical harmonics [21] 
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' 21 + 1 

Y i : n (n b )Y lm (n a ) = ^—P l (n a -n b ) (159) 

m=— Z 

and where Pi is a Legendre polynomial: 

A. 2 System response algebra 

From equation ( |103|) , i.e., 
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^afe + r/ E 
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(s 2 + 2)( s 2 + ^ ) 
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(nl) 
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s 2 + n 2 
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(161) 

we must first isolate cjb(s), then find inverse Laplace transforms to revert to time domain quantities. 
Substituting the values of w (s) and £® w (s) from fll54.a|) and (Pi ) into (|l~6ll) we find 



q a (s)= Y, ^ m) (s)g (lm) (s) 



a = 1 J 



(162) 



Z = and 2 
m = — L 



where 



s 2 + ft 2 
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<5 a 6 + E 
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ft 2 S 



( s 2 +0 2) (s 2 + ^) 



(nl) 
2 -^afe 



^m(nfe) 

(163) 

Now, using the convolution theorem of Laplace transforms, we see that the time domain version of 
equation (|162| ) is 



= E f $ i M (* - <7 (M (0 ^' , a = 1, . . . , J 



(164) 



where $i' m ^(t) is the inverse Laplace transform of ( |163| ). The inverse Laplace transform of &a\s) can 
be expediently calculated by the residue theorem through the formula [26| 



<jj(Zm)(£) _ 2-xi ^ {residues of ^ m \s) e st at its poles in complex s - plane j 



(165) 



Clearly thus, the poles of &a (s) must be determined in the first place. It is immediately clear 
from equation ( |163j ) that there are no poles at either s = 0, or s = ±i$7, or s = ±iuj n i, for there are 
exactly compensated infinities at these locations. The only possible poles lie at those values of s for 
which the matrix in square brackets in (|163|) is not invertible, and these of course correspond to the 



zeroes of its determinant, i.e., 



A(s) = det 



6ab + 71 ^Tw 



K ab {s) 







poles 



(166) 



There are infinitely many roots for equation ( |166| ), but analytic expressions cannot be found for 
them. Perturbative approximations in terms of the small parameter n will thus be applied instead. It 
is assumed that 
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Figure 7: Position of the roots ( [169D in the complex s-plane for SI tuned to the lowest frequency of 
the sphere's spectrum. Note that s n i is closer to iuj n i than sq± are to iSl; this is a consequence of the 
differences being proportional to n and n 1 ^ 2 , respectively, as seen in equations (|169|). 



SI = u no i (167) 
for a fixed multipole harmonic {nolo}. Equation ( |166| ) can then be recast in the more convenient form 
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(168) 



Since 77 is a small parameter, the denominators of the fractions in the different terms in square 
brackets in ( |168| ) must be quantities of order r\ at the root locations for the determinant to vanish at 
them. A distinction however arises depending on whether s 2 is close to — ft 2 or to the other — u; 2 z . 
There are accordingly two categories of roots, more precisely: 



si 



„2 
s nl 



Si 2 (l + X 1 r/ 1/2 + Xi + . 
( 1 + b\ n + b\ ' t] + ... 
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(169.a) 
(169.b) 

(nl) 



Figure ^ provides a visual characterisation of these roots. The coefficients xi> Xi>— an d ^1 

&i can be calculated recursively, starting form the first, by substitution of the corresponding series 
expansions into equation (|168| ). The lowest order terms are easily seen to be given by 
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(170) 
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and 



det 



9 2 



u >nl i>0 A Anl) 
i O l Oab - X a b 



nl 



(171) 



respectively. Both equations (170) and ( |171| ) are algebraic eigenvalue equations. As shown in ap- 
pendix [A.3| , the matrix xi& nas & t most (21 + 1) non-null positive eigenvalues — all the rest up to J 
are identically zero. 

As a final step we must evaluate (|165|) . This is accomplished by standard textbook techniques (see 
e.g. [^]); the algebra is quite straightforward but rather lengthy, and I shall not delve into its details 
here, but quote only the most interesting results. It appears that the dominant contribution to <l>a (i) 



comes from the poles at the locations ( |169.a ), whereas all other poles only contribute as higher order 
corrections; generically, $a\t) is seen to have the form 



$« m \t) oc rT 1/2 (sinu; c+ t - sinu; c _t) S llo + O(0) 
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In Laplace domain one has, 
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(174) 



Detailed calculation of the residues [57] yield equation ( |106| ), which must be evaluated for each 
particular tuning and resonator distribution, as described in section ||. 



A. 3 Eigenvalue properties 

This Appendix presents a few important properties of the matrix P/(n a -nt) for arbitrary / and resonator 
locations n a (a=l,. . . ,J) which are useful for detailed system resonance characterisation. 



Recall the summation formula for spherical harmonics [21 
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where Pi is a Legendre polynomial 
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2 1 n dz l 

To ease the notation I shall use the symbol V\ to mean the entire JxJ matrix P;(n a -nb), and 



introduce Dirac kets \m) for the column J- vectors 

*lm( n l) 



\m) 



I 4vr 
21 + 1 
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(177) 



These kets are not normalised; in terms of them equation ( |175| ) can be rewritten in the more compact 
form 
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I 

Vi = \ m )( m \ ( 178 ) 

m=— I 

Equation ( |178| ) indicates that the rank of the matrix V\ cannot exceed (21 + 1), as there are only 
(21 + 1) kets \m). So, if J > (21 + 1) then it has at least (J — 21 — 1) identically null eigenvalues — there 
can be more if some of the n a 's are parallel, as this causes rows (or columns) of Vi to be repeated. 

We now prove that the non-null eigenvalues of V\ are positive. Clearly, a regular eigenvector, \<j)), 
say, of Vi will be a linear combination of the kets |m): 

I 

Vi \4>) = C 2 !</>> , \<t>) = £ <Pm \m) 

m=—l 

where Q 2 is the corresponding eigenvalue, having a positive value, as we now prove 
is substituted into the first then it is immediately seen that 

I 

£ (C 2 <W - Hm')) <Pm> = (180) 

m'=-l 

which admits non-trivial solutions if and only if 

det (C 2 <W - {m\m')) =0 (181) 

In other words, Q 2 are the eigenvalues of the (21 + l)x(2Z + 1) matrix (m\m'), which is positive 
definite because so is the "scalar product" (0| </>'). All of them are therefore strictly positive. 
Finally, since the trace is an invariant property of a matrix, and 

J J 
traced) = £ P/(n a -n a ) = £ 1 = J (182) 

a=l a=l 

we see that the eigenvalues ( 2 add up to J: 

J 

traced) = £ £ = £ (a = J (183) 
a=l Ca^O 
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